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INVESTIGATION OF NORTHROP F-5A WING BUFFET INTENSITY 
IN TRANSONIC FLIGHT 


By Chintsun Hwang and W.S. Pi 

Northrop Corporation, Aircraft Division 
Hawthorne, California 


SUMMARY 


The report describes the findings and conclusions of a buffet flight test and data 
processing program entitled "Investigation of Buffet Intensity in Transonic Flight." 
The program used a fully instrumented Northrop F-5A aircraft to acquire buffet 
pressure and response data during transonic maneuvers. In the report, the detailed 
instrumentation, the data flow sequences, and the flight conditions used in the maneu- 
vers are described. The dynamic buffet pressure and response data, as recorded by 
Pulse-Code Modulation (PCM) and FM systems, are presented in real-time format for 
general evaluation. These are followed by spectral and statistical processing pro- 
cedures where the pressure and response spectral density functions and other 
statistical data are presented. Furthermore, the spatial correlation of the individual 
pressure data and the correlation between buffet pressure and aircraft response are 
described. Also covered is an analytical procedure to correlate the aircraft response 
data with computed responses based on measured buffet pressures. 
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Section 1 


INTRODUCTION 

Buffet is a dynamic behavior of aircraft flying at a high angle of attack. It 
occurs in the transonic region where the aircraft has substantial forward speed, and 
the angle of attack becomes large due to certain maneuvers initiated by the pilot. The 
initial phase of buffet is called buffet onset when noticeable vibrations or oscillations 
are observed. When the angle of attack is further increased, the aircraft responds to 
the dynamic buffet pressure with severe structural vibration which may be coupled with 
rigid body motion oscillations or control system -induced oscillations. This latter 
phenomenon of sustained, large-amplitude aircraft buffeting is called wing rock. 


1 . 1 Background 

As is well known, the major cause of aircraft buffet in the transonic region is 
flow separation on top of the wing surface when the angle of attack reaches a certain 
amplitude. The separated flow may be induced by a shock or by other disturbances. 
The dynamic loads exerted on the aircraft may be due to the instability of the shock 
(in intensity and location), or due to the dynamic pressure components of the turbulent 
flow especially in the separated region, or both. Furthermore, the separated flow 
wake may engulf the tail sections, which causes additional instability in the aircraft 
responses and deterioration of the flight control qualities. 

Earlier work on aircraft buffeting may be traced to the associated subject of 
stalling flutter (References 1, 2, 3). Both stalling flutter and aircraft buffeting are 
caused by flow separation at high angle of attack. In stalling flutter, the airfoil motion 
(for instance, torsional motion) contributes to the cyclic flow pattern which causes a 
sustained oscillation. In aircraft buffeting, the structural elasticity usually is not a 
major contributing factor (at least in the initial phase of the transonic maneuver). The 
buffeting loads are aerodynamic in origin; they exist even if the aircraft, or aircraft 
model, is perfectly rigid. 
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The subject of aircraft buffeting started to attract the attention of designers and 
research workers in the 1950's when aircraft development reached such a stage that 
transonic flights and maneuvers became a routine matter. Typical works by Huston 
and associates appearing at that time are given as References 4 through 7. Wind 
tunnel investigations on airfoil or aircraft buffet in this period include those by 
Humphreys, Coe, Sutter, and associates (References 8, 9, 10, 11). A corresponding 
work carried out in England by Pearcey is given as Reference 12. Pearcey at a later 
date also wrote a very informative article on shock -induced separation and buffet 
loads (Reference 13). 

Starting in 1960, research and development on aircraft buffeting were conducted 
at an accelerated pace. A large number of papers and reports were published dealing 
with the various aspects of aircraft buffet including theoretical treatment (e.g. , 
References 14, 15), flight test evaluation (References 16, 17, 18), and wind tunnel 
testing (References 19 through 22). Still more recent works dealing with buffet of 
present-day and future aircraft are listed as References 23 through 30. Reference 26 
describes the results of F-5A buffet flight tests, with emphasis on stability and con- 
trol. Comprehensive buffet onset and wing tuft data and limited amounts of total and 
static pressure data were acquired in this test program. Reference 27 gives the 
buffet onset data of the F-8D based on both flight tests and wind tunnel tests. Refer- 
ences 28 through 30 describe the buffet characteristics and the dynamic buffet loads 
of the YF-4E aircraft where the dynamic pressure components on the top surfaces of 
the wing were greatly affected by vortices generated at the leading edges and wingtips. 


1.2 Basic Program Approach 

The work described in this report deals with the flight test and data evaluations 
of the Northrop F-5A aircraft in transonic maneuver. A unique feature of the subject 
program involves the extensive dynamic pressure instrumentation using miniaturized 
semi-conductor type transducers. Both fluctuating buffet pressure and significant air- 
craft response data were recorded on a wideband FM recording system used for 
real-time and spectral analyses. For instance, using digital tapes generated from 
the precisely clocked FM data, spatial correlations (as a function of the time 
parameter) and normalized cross-covariance functions were generated between two 
pressure stations, between the pressure and the response function, and between two 
selected response functions. 
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Another feature of the subject program involved the spectral processing of the 
flight data using various time segments obtained during the transonic maneuver. The 
basic transonic maneuver was the wind-up turn executed at specific Mach number and 
starting altitude. The complete maneuver lasted from 12 to 18 seconds which included 
level flight entry, buffet onset, wing rock, recovery, and the intermediate stages. 

The data acquired in a small time segment of the maneuver (approximately two 
seconds, depending on frequency range) were assumed to be stationary so that auto- 
and cross-power spectra could be generated. The display of a set of power spectra 
generated for various time segments of the maneuver yielded results that illustrated 
the various phases of the buffet phenomenon. 

To complement the flight test program, work was carried out to correlate the 
response data using the measured buffet pressure and analytical computation tech- 
nique taking into account the aircraft structural and aerodynamic characteristics. 

The correlation procedure and other statistical processing techniques are described 
in a later section of the report. 

In summary, the data evaluation and correlation work of the subject program 
consisted of the following specific tasks: 

1. Make cross-spectral analysis of the wing pressures and compute 
the wing responses. 

2. Make spectral analysis of the response parameters for comparison with 
the computations and with each other. 

3. Determine the probability distribution of one of the wing response para- 
meters. 

4. Make other analysis as may appear useful in estimating buffet loads 
for general use. 

The program was conducted under the sponsorship of Ames Research Center, 
National Aeronautics and Space Administration. The NASA program monitor was 
Mr. Donald Buell. At Northrop, Messrs. Carl Balmer and Merv Burne assumed 
major responsibilities in aircraft instrumentation and data recording. Mr. Frank 
Jacobs contributed to the digital data processing. Dr. Wilford Wong carried out the 
evaluation of the real-time data and contributed to the preparation of the report. The 
authors wish to acknowledge the assistance and support rendered by Northrop F-5 
project personnel inbdfn Dynamics Department and Aerodynamics Department. The 

F-5A aircraft bailment was arranged through Major Horace Russell of the Air Force 
Flight Dynamics Laboratory. The authors also wish to' acknowledge the enthusiastic 
support of Mr. Charles Coe of N/lSA and Mr. Hugo Pink of Northrop. 
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1. 3 Symbols and Units 

Fourier transform of the deterministic time function 
e^ y seconds 

Span (clean tips), meters 

Reference semi-chord, meters 

Equivalent resolution bandwidth, Hz 

Mean aerodynamic chord, meters 

Aerodynamic derivatives as defined in Chapter 6, 
Reference 59 


Deterministic time functions for various time segments 
Frequency, Hz 
Cutoff frequency, Hz 

Fourier transforms, newton/m eter2, etc. 

Modal structural damping coefficients 

(1-/9) • (Modal aero damping coefficients) + 1 

2 

Gravitational acceleration, meters/ second 
Flight altitude, meters 

Modal influence functions for various time segments, 
G/newton, etc. 

Reduced frequency 
Maximum correlation lag, seconds 
Total mass of the aircraft, kilograms 
Modal masses, kilograms 
Mach number 

Total number of time segments in Appendix III 
Sample size 

Longitudinal, lateral, and normal accelerations, multiples 
of G 



p 


p 

2P 

2PX 

SPY 

Q < i 

r 



t, T 


)a 

)b 


AT ~\ )b" \ )a 
V 


V 

c 

w 


x, y, z 


a 

P 


P, 


) 


y 


2 

( ) 


2 

Probability density function, meter /newton, etc.} also 
pressure load, newtons/meter 

Cumulative probability distribution function 

Total force, newton 

Total torsional moment, newton-meter 

Total bending moment, newton-meter 
Segmentwise stationary modal forces, newtons 

Distance between two points, meters 

Correlation functions, (newton/m eter^) , etc. 

Reference length in Appendix I, meters; also reference 

2 

area, in Appendix III, meter 
Time, seconds 

Initial time of various time segments, seconds 
Terminal time of various time segments, seconds 
Time span, seconds 

Time span for each individual time segment, seconds 

Aircraft speed, meters/second 

Convection speed, meters/ second 

Structure deflection at a specified location, meters 

Cartesian coordinates, meters 

Truncated function of x(t) with time span T 

Angle of attack, degrees 

A faetor defining the portion of induced pressures to be 
assigned to the upper wing surface 

(1- P ) • (Imaginary part of modal aero-damping 
coefficients) 

Coherence functions 
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Spatial decay constant 



e 

€ 

( ) 
t 

"( ) 

6 

P 

p o 

P ( ) 

*( ), 0 ( ) 

0 -, , 

'< ) 
w= 2ni 

"( ) 

Matrix Conventions 

40 

;|fcj 

- L J 
[ ] T 
[ I'" 


Leading edge flap angle, degrees 
Trailing edge flap angle, degrees 
Normalized standard error 

(1-/?). (Real part of modal aero-damping coefficients) 

Relative magnitude of the induced pressure at the top wing 
surface to the measured pressure 

Mean values 

Phase angle, radians 

3 

Air density, kilograms/meter 

Normalized auto- correlation function 

Correlation function coefficient or normalized 
cross-covariance functions 

2 2 

One-sided spectral density functions, (newton/meter ) /Hz 
Raw estimate of the true power spectral density function 

£) O 

0^ y (newton/meter ) / Hz, etc. 
rms values 

Circular frequency, radians/ second 

Non-zero natural circular frequency, radians/second 


Square or rectangular matrix 
Column matrix 
Diagonal matrix 
Row matrix 
Transposed matrix 
Inverse matrix 
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Matrices 


[A] 

[ A F ] 

?<: 

[c] 


rc F j 


r%j 


[H] 

r. ( ,i 

ri n J 

r* F J 

[ k ( ,] 
Cm f J 

[M () ] 

[«,,] 


Wing subarea matrix associated with pressure 
transducers, meter 2 

Generalized AIC matrix, kilograms 

Flutter characteristic matrix, newtons/meter 

Generalized work coefficient matrices due to external 
forces 

Damping coefficient matrix corresponding to control 
point velocity, newton-second/ meter 

Oscillatory AIC matrix for control point motion 

Generalized damping matrix, newton-second/meter 

Fourier transforms of modal force matrix, newtons; also 
Fourier transforms of collocation force matrix in 
Appendix I 

Structural damping coefficient matrix 

Fourier transforms of deflection matrices corresponding 
to control points or reference point, meters or radians 

Modal transfer function matrix, meters/newton; 
or radians/ (newton-meter)- 

Modal deflection influence function matrix for various 
time segments, meters/newton or radians/newton-meter 

Modal acceleration influence function matrix for various 
time segments, newton \ etc. 

Generalized stiffness matrix, newtons/meter 

Structural stiffness matrices, newtons/meter 

Generalized (modal) mass matrix, kilograms 

Mass matrix, kilograms, etc. 

Generalized mass plus aero matrices, kilograms, etc. 

Measured pressure matrix for various time segments, 
newtons/ meter 2 
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[Q () ] 

[s, )] 

[t] 


[w] 

[X] 


[Q f ] [Q r ] 


Modal shape matrices corresponding to pressure 
transducer locations 


Two-sided spectral density matrices, 
(radian/second), etc. 


(newton/ meter^ V 


Rigid modal shape matrix corresponding to control points, 
meters 

Aeordynamic weighting matrix 

Modal shape matrix corresponding to pressure transducer 
locations, non-dimensional or meters 


c x <>] 

L Yj 

0 ]-[ H ]- 1 



Flexible modal matrix correspondings to control points or 
reference point, non-dimensional or meters 

Modal shape matrix corresponding to a specified accelerom- 
eter location, non-dimensional or meters 

Modal impedance matrix, newtons/meter or newton-meters/ 
radian 

Impedance matrix corresponding to the downward displace- 
ment and the pitch angle, newtons/meter, newton-meters/ 
radian, etc. 

Fourier transforms of modal amplitudes (generalized 
coordinates), meters or radians 

Error in series solution 

Eigenvector matrix corresponding to control points and 
reference point, meters or radians 


Subscripts and Signs 


( ) 
( ) 
( ) 
( ) 

( ) 


1 

2 

b 

f 

F 


Refers to control points (Appendices I and HI) 

Refers to reference point (Appendices I and III) 

Refers to buffet pressures 

Refers to modal forces; also refers to flexible modes in 
Appendix III 


Refers to flexible modes 



< >i 

Refers to induced aerodynamic pressures 

< >k, ( >1 

Refers to modes in Appendix III 

< >m 

Refers to measured pressures 

( >P 

Refers to pressures 

( ) r ,( ) s 

Refers to time segments in Appendix III 


Refers to rigid body modes 

( .)« 

Refers to modal amplitudes 

1 1 

Denotes absolute value 

(.*) 

Denotes time derivative 

* 

( ) 

Denotes complex conjugate 

(~ ) 

Denotes complex mass which includes oscillatory 


aerodynamic effects 

T ) 2 

Denotes mean square value 

A( ) 

Denotes increment 

Special Symbols 


(N/M**2)**2/HZ 

(newton/ meter 2 ) 2 / Hz 

RW32B18 

Strain gage at RHW.S. 32 designated B18, used with other 
gages to determine sectional bending moment 

RW32B24 

Strain gageat RHW.S. 32 designated B24, used with other 
gages to determine sectional bending moment 

RW32S21S 

Strain gage at RH W.S. 32 designated S21S, used with 
another gage to determine sectional shear force 

RW32S44S 

Strain gage at RH W.S. 32 designated S44S, used with 
another gage to determine sectional shear force 

RW68B18 

Strain gage at RH W.S. 68 designated B18, used with other 
gages to determine sectional bending moment 

RW68T30 

Strain gage at RH W.S. 68 designated T30, used with other 
gages to determine sectional torsional moment 
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Section 2 


FLIGHT TEST PROGRAM 

The buffet flight test program of the F-5A was carried out on a fully 
instrumented aircraft (No. 63-8372, N6009). This is the same aircraft used in a 
previous buffet test program which collected excellent data on buffet onset, wing sur- 
face tufts behavior, and other items of general interest (see Reference 26). As men- 
tioned in Section 1, the objectives of the present program include the acquisition and 
processing of the dynamic pressure and response data covering all stages of the trans- 
onic maneuver. For this purpose, new and additional instrumentation was installed 
on the aircraft prior to the flight tests. Because the buffet pressures in the shock and 
separated flow region are at times quite transient in nature, and because the pressure 
spectrum covers a wide frequency range corresponding to all stages of the transonic 
maneuver , an FM recording system was installed on the aircraft to record the physi- 
cal phenomenon. 


2.1 Aircraft Configuration 

The F-5A used in the test program was a single-seat fighter capable of carrying 
stores at wing fuselage pylon stations. Due to the requirements of another research 
program being carried out simultaneously on the same airplane, the test was conducted 
with two wingtip stores (AIM-9B Missiles) with guide rails; otherwise the wing was 
clean. A three-view drawing of the F-5A is shown as Figure 1. 

The basic wing airfoil of the F-5A is a NACA 65A004.8 section modified as 
follows. From the 40% chord to the 66.6% chord, the thickness distribution corre- 
sponds to a NACA 004.8-64 airfoil. Aft of the 82% chord, the airfoil is composed of 
straight lines forming a 7° trailing edge angle. Between the 66.6% chord and the 82% 
chord, the section has been curve-fitted so that the ordinates, slopes, and curvature 
at the end points are continuous. The airfoil is cambered to a NACA (. 65) 50 five- 
digit series mean line. The mean line ordinates are obtained by multiplying the stan- 
dard NACA 250 series mean line by the ratio .65/2. The mean line is rotated so that 
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the straight aft portion lies in the wing chord plane and the camber appears as a 1% 
leading edge droop. The airfoil is constant over the entire span. 


For the flight test program, the leading edge flaps were fixed at 0° for some 
flights and at 4° for others. Likewise, the trailing edge flap systems were modified to 
provide flap travel limits of 0° to 8° for certain flights and 0° to 12° for others. (The 
flap could not be stopped at intermediate positions.) The basic wing airfoil section 
with variable flap configurations is shown in Figure 2. 


Data Description of Wing 


Area (Reference) 

15. 79 m 

(170.00 ft ) 

Span (clean tips) 

7.696 m 

( 25.25 ft ) 

Aspect Ratio 

3.75 


Taper Ratio 

.20 


Sweepback (25% Chord) 

24° 


Mean Aerodynamic Chord (Reference) 

2. 356 m 

( 7. 73 ft ) 

Thickness Ratio 

4.8% 


Dihedral Angle 

0 


Incidence Angle 

0 


Root Chord 

3.420 m 

( 11.22 ft ) 

Tip Chord 

. 683 m 

( 2.24 ft ) 

15% Chord Leading Edge Flap (with 
built-in twist) Deflection When 
Wholly Extended (Nominal flap 
deflection angle is measured at 
the root) 

23° at the root 
18° at the tip 


T.E. Flap Deflection When Wholly 
Extended 

20° 


AIM-9B Missile Diameter 

0.127 m 

( 5.00 in ) 

AIM-9B Missile Length 

2. 832 m 

( 9.29 ft ) 


The horizontal tail of the F-5A is all-movable and has an unmodified NACA 
65A004 section with the following physical dimensions and configurations: 


Data Description of Horizontal Tail 

Area, Total 
Area, Exposed 
Aspect Ratio, Exposed 
Sweepback (25% Chord) 


5.48 m 2 ( 59 ft 2 ) 

3.07 m 2 ( 33.03 ft 2 ) 

2.88 
25° 
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Taper Ratio, Exposed 
Cathedral Angle 

Mean Aerodynamic Chord, Exposed 


.33 

4° 

1.122 m 


( 3.68 ft) 


The vertical tail airfoil is a modified NACA 65A004 section. The airfoil is com- 
posed of straight lines aft of the 72% chord. Between the 60% chord and 72% chord, 
the airfoil consists of a transition curve which matches the ordinate, slope, and 
radius of curvature of the true airfoil at 60% and the ordinate, slope, and infinite 
radius of curvature at the 72% chord. 


Data Description of Vertical Tail 


Area, Exposed 

3.83 m 

Aspect Ratio, Exposed 

1.22 

Taper Ratio, Exposed 

.25 

Sweepback (25% Chord) 

25° 

Mean Aerodynamic Chord, Exposed 

1.990 m 


Useable Fuel, Weight, and Power Plant Data 


(41.2 ft 2 ) 


( 6. 53 ft ) 


Internal Fuel 1, 719 KG, 2207 liters (3, 790 lbs, 583 gal.) 

Aircraft Weight with Full 

Fuselage Tank and AIM-9B's 6, 211 KG (13, 694 lbs) 

AIM-9B Missile Weight 

(Including Launch Rail), Each 92.5 KG ( 204 lbs) 

Power Plant 2 (J85-GE-5) Turbo-jets with Afterburners 


2. 2 Test Aircraft Instrumentation 


In the flight test program, acquisition of data, other than flight condition param- 
eters, was accomplished by a magnetic tape system utilizing digital pulse code modu- 
lation (PCM) and analog frequency modulation (FM) recording. . The PCM data were 
recorded in serial format on two tracks of a 14-track, 2.54-cm-wide (1-inch) mag- 
netic tape. Eleven tracks were utilized for FM recording and one track was used to 
record pilot's voice annotations and event mark signals. The wideband magnetic tape 
recorder, operating at 38.1 cm per second (15 inches per second), provided suitable 
frequency response characteristics and allowed 60 minutes of recording time on one 
tape reel. The FM data were subsequently digitized at a sample rate of 5,000 points 
per second to be used for real-time and spectral processing. 
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The PCM instrumentation installed in the test aircraft consisted of two encoding 
units, analog signal conditioning and presample data filtering equipment, synchro-to- 
digital and frequency -to-digital converters, a digital control unit, and a subsystem 
interface unit. The PCM encoders was used to perform analog signal multiplexing, 
analog-to-digital conversion, and digital output formatting. The system operated at 
a rate of 100,000 bits per second to output successive frames of 72 10-bit digital 
words (± 512 counts full-scale) in serial format at a rate of 138.9 frames per second. 
Direct entry of digital words from the synchro and frequency converters was 
accomplished by the digital control unit. The nominal frequency response of the PCM 
system operating at 138.9 samples per second was 30 Hz as established by the pre- 
sample data filters which rolled off at 18db per octave above 30 Hz to prevent aliasing 
of encoded data. Five parameters were recorded at a rate of 555.6 samples per 
second by cross-strapping together four prime channels each. These channels were 
provided with 120-Hz filters. All parameters recorded by the PCM system specifi- 
cally for the buffet flight test program are listed in Table 1 . 

Since the buffet test program was carried out simultaneously with another flight 
test program, additional PCM data were recorded which were available to the present 
program. This point will be explained later relating to flight test data flow. 


17 



TABLE 1 


PULSE CODE MODULATION RECORDING 


IH 





Samples 

Item 


Parameter 


Per 

IBM 





Second 

■ 

Pressure Pickup No. 4, Upper RH W.S. 128.31,* 90%C (85% b/2) 

556 

2 


19, Lower RH W.S. : 

1.28.31 , 20%C (85% b/2) 



3 


20, 


60%C 




4 


21, 


90%C 




5 


24, 

W.S. 

93.00, 60%C (61% 

b/2) 



6 

Normal Load Factor at C.G. 




1 

39 

7 

Longitudinal Load Factor at C.G. 





8 

Normal Load Factor, RH Wingtip, Aft 

80%C 




9 


RH Wingtip, Fwd 

30%C 




10 


Pilot's Seat 





11 

Lateral Load Factor, Pilot's Seat 





12 

Angle of Attack, CADC (Central Air Data Computer) Output 



13 

Pitch Rate 







14 

Pitch Attitude 






15 

Right-Hand Aileron Position 






16 

Strain Gage, RH W.S. 127, Designated B40 





17 

Right-Hand Trailing Edge Flap Hinge Moment 


1 


18 

Camera Frame 






LiL 

Event Mark 






* Multiply Wing Station data by .0254 to obtain Station distance in meters. 


An FM instrumentation system permitted recording of dynamic data at frequen- 
cies up to 1, 600 Hz based on a modulation index of 5. Frequency multiplexing was 
used to carry three data parameters on each of the eleven FM recording tracks. The 
center frequencies of the discriminators were 64, 96, 128 KHz respectively. The 
bandwidth of the modulated frequency was ± 8 KHz. The parameters recorded in this 
manner are listed in Table 2. The details of the pressure transducers (PCM items 
1 through 5, FM items 1 through 20) are described in a latter part of this subsection. 
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TABLE 2 


FREQUENCY MODULATION RECORDING 


Item 


Parameter 


21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 


1 

Pressure 

Pickup No. 1. Upper RH W.S. 

128.31,* 

20%C 

(85% 

b/2) 

2 


2, 



40%C 



3 


3, 



60%C 



4 


4, 



90%C 



5 


5, 

W.S. 

110.30. 

20%C 

(72.8%b/2i 

6 


6, 



40%C 



7 


7, 



60%C 



8 


8, 



75%C 



9 


9, 



90%C 



10 


10, 

W.S. 

93.00, 

20%C 

(61% 

-b/2) 

11 


11, 



40%C 



12 


12, 



60%C 



13 


13, 

W.S. 

71.40, 

60%C 

(47% 

b/2) 

14 


14, 

W.S. 

71.40, 

80%C 

(47% 

b/2) 

15 


15, 

W.S. 

50.00, 

20%C 

(33% 

b/2) 

16 


16, 



40%C 



17 


17, 



60%C 



18 


18, 



80%C 



19 


22, 

W.S. 

L26.31, 

75%C 

(83. 4% b/2) 

20 


23, 

W.S. 

71.40, 

40%C 

(47% 

b/2) 


Strain Gage , RH W.S. 32 , 


Designated B24 
Designated S21S 
Designated S44S 
Designated B18 
RH W.S. 68, Designated B 18 
RH W.S. 68, Designated T40 
Normal Load Factor , RHWingtip, FwdW.S. 144.85 

RHWingtip, Aft W.S. 144.85 
C.G. 

LH Wingtip, Aft 
Right-Hand Aileron Hinge Moment 
Event Mark 


30%C 

86%C 


Multiply Wing Station data by . 0254 to obtain Station distance in meters. 


Flight and engine condition instruments and aircraft system status lights were 
located on a panel and photographed by a 35mm frame-pulsed camera at one-second 
intervals. The parameters recorded are listed in Table 3. 
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TABLE 3 


PHOTO-PANEL RECORDING 


Item 

Parameter 

1 

Indicated Airspeed, Production Boom 

2 

Pressure Altitude, Production Boom 

3 

Outside Air Temperature, Rosemount Probe 

4 

Left-Hand Fuel Quantity 

5 

Right-Hand Fuel Quantity 

6 

Event Light 

7 

Correlation Counter 


2.3 Dynamic Pressure Transducers 

To record the dynamic pressures, 24 miniaturized differential pressure trans- 
ducers (SENSOTEC Model SD-M-7F) were installed on the right wing of the test air- 
craft. These pressure transducers had a diaphram sensor unit. Four-arm semicon- 
ductor gages connected in a bridge circuit were mounted on the internal face of the 
diaphram. Temperature compensation was achieved through thermistors and active 
temperature components (resistors) located within the transducers. 

In addition to the four-wire bridge connection , each transducer had a reference 

pressure outlet which was connected to the reference source. The transducers 

selected for the test program had a diaphram diameter of 0.635 cm (0.250 inch) with 

0.076 cm (0.030 inch) thickness and 1.270 cm (0.500 inch) overall length. The 

2 

nominal pressure range was 1. 724 N/cm (2.5 psi) with an overpressure tolerance of 
2 

3.448 N/cm'". In addition to the calibration data furnished by the manufacturer, the 
transducers were calibrated dynamically prior to installation using Photocon Acoustic 
Calibrator and B & K Microphone. 

To minimize the interference to natural flow conditions, the transducers were 
flush-mounted on the right wing surfaces together with a 0. 1016 cm (0. 040 inch) plexi- 
glass jacket assembled in 30.48 cm (12 inch) strips. The plexiglass jacket covered 
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almost 80% of the exposed area of the top surface and was laid in two 30.48 cm (12 inch) 
strips on the bottom surface. (The covered areas are indicated by the heavy outlines 
on Figure 3.) Prior to the jacket installation, all the transducer harness wires and 
reference pressure tubings were assembled and laid on the wing surface, and the ter- 
minals were led into the internal space of the wing structure through rivet holes. 
Properly located channels were cut on the underside of the plexiglass jacket, together 
with circular holes, prior to its bonding to the metal skin of the wing. These channels 
served to accommodate the transducer harnesses, while the circular holes were lo- 
cated where the pressure transducers were to be mounted. A 15-to-l taper was cut on 
the outer circumference of the jacket to minimize flow interference. After the plexi- 
glass sheets were trimmed and bonded by adhesives to the wing, the transducers were 
cemented in place and connected to the electric wirings and the reference pressure 
tubing (Figure 4). Pre-cut metal plugs were cemented over the remaining space of the 
circular hole to form a smooth airfoil surface on the wing. A photograph showing the 
top surface of the right wing after the installation of the pressure transducers is 
shown in Figure 5. 

Internally, the reference pressure tubings from all transducers were plumbed 
into a common manifold chamber located inside the lower right hand section of the 
fuselage at F.S. 349.5*, where a precision altimeter also was installed. The pressure 
transducers, numbered 1-18, 22, and 23 on the top wing surface (see Figures 3 and 5), 
were listed as FM recordingdata entries in Table 2. The four pressure transducers 
at the bottom surface, designated Numbers 19, 20, 21, 24, were shown in Table 1 as 
PCM data entries. Pressure transducer Number 4 was connected to both FM and PCM 
systems to evaluate the recording system accuracies. 

The signals from the pressure transducers were amplified through a two-stage 
amplifier unit with adjustable gains (BURR-BROWN Models 3620K-3621K). The signal 
conditioning circuit included a high-pass filter with l/2-Hz cut-off frequency between 
the two amplifier stages. The signal from the first stage was input to the PCM system 
through a presampling (low-pass) filter. The signal from the second stage was input 
to the FM system so that the low-frequency and steady-state pressure data were not 
reflected in the FM data. A schematic showing the FM recording system for the upper 
Wing surface pressure transducers is given in Figure 6. 


♦Multiply Fuselage Station Number by .0254 to obtain Station distance in meters. 
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UPPER 

SURFACE 


LOWER 

SURFACE 


On F-5A Right Wing 
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23 


Figure 4. Sensotec Pressure Transducer Mounted on the F-5A Wing 




Figure 6. Recording Circuit-Upper Wing Differential 
Pressure Transducers 
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2.4 Data Flow and Processing Sequence 


As mentioned previously, both FM and PCM recording systems were used in the 
flight test program . Figure 7 shows the PCM and telemetery data flow from test air- 
craft through the Edwards Air Force Base Test Facility to the Northrop Hawthorne Data 
Processing Facility. Except for differences in signal conditioning amplification and 
multiplexing techniques, as well as the final decoding and digital tape conversion, the 
FM data were handled in a similar flow pattern as the PCM data. Certain flight data 
parameters were recorded in addition to those listed in Table 1 (see Figure 7). These 
additional data were shared with a simultaneous flight test program and were used in 
the final buffet data evaluation. As part of the data flow concept, the Edwards Data 
Processing System provided quicklook data editing and instrumentation diagnostics 
from the flight data tape. 

Upon completion of a test flight, the PCM data were played back from the air- 
borne data tape. The recorded high-frequency analog channels were read out into FM 
discriminators which separated out each VCO frequency and demodulated the analog 
data. These high-frequency analog data were then recorded on oscillographs for pre- 
liminary editing operations. Further analysis of the analog data was performed on 
selected data segments at the Hawthorne Data Processing Facility. The Hawthorne 
Facility was equipped with tape playback machines, discriminators, digital-to-analog 
converters, and digitizing and transcribing equipment to process the buffet test data. 
Quick evaluation of PCM and FM recording was performed by transcribing the data on 
an oscillograph tape to determine periods of interest for further processing. The data 
flow to the computer and to the SC-4020 system is shown in Figures 8 and 9, 
respectively. Photo recorder data which were eventually integrated into the other 
recorded data flow also were reduced. 
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Figure 7. Flight Test Data Flow Diagram 
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FIGURE 8. DATA FLOW TO COMPUTER FOR BUFFET TEST PROGRAM 



Figure 9. Data Flow to SC 4020 for Buffet Test Program 
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2.5 Flight Conditions 


The basic maneuver used in the test program was the wind-up turn at constant 
Mach number. Specifically, wind-up turn maneuvers to maximum lift were performed 
corresponding to the configurations and flight conditions of Table 4. Prior to each 
maneuver, the aircraft was held in level flight at the desired test condition for 60 
seconds or more to stabilize the transducer reference pressure, etc. Data were 
recorded from 1-G level flight through recovery from the maneuver. The speed brakes 
remained closed and the wing flap settings were maintained throughout each maneuver. 
The tests were conducted in a manner such that full aft stick was attained for 7 to 12 
seconds after maneuver initiation. Satisfactory data were obtained from at least two 
maneuvers for each condition. Furthermore, roller-coaster maneuver from -0.5G to 
1. 8G load factor was performed during each flight at 7,620 m (25, 000 ft), and 0. 70 
Mach number to provide the zero-G strain gage response references necessary for 
computation of absolute values of W.S. 32 bending moment, shear, and torsion loads. 

The altitude-Maeh number combination (Table 4) chosen for the wind-up turns 

3 2 

was selected based on a constant q of approximately 14. 36 x 10 N/M (300 psf). The 
three conditions are plotted on the speed profile of the F-5A corresponding to various 
gross weight figures (see Figure 10). 


TABLE 4 

TEST CONFIGURATIONS AND FLIGHT CONDITIONS 
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Figure 10. Speed Profile of F-5A with Maximum Thrust and (2) AIM-9B 
Missiles (Buffet Test Conditions Are Marked by Triangles) 
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SECTION 3 


BUFFET DATA ACQUISITION AND INITIAL PROCESSING 

Following extensive ground checkout of test aircraft new instrumentation, includ- 
ing the dynamic pressure transducers, the associated signal conditioning and amplifier 
circuits, the FM recording system, etc., shakedown flights were made in late 1971. 

The actual buffet flight tests were conducted starting in February 1972. All the flight 
conditions specified in Table 4 were performed and the data recorded. After a fast 
scan of flight test results, the recorded data were converted into digital tapes which 
were subjected to both real-time and spectral processing. In this section, the process- 
ing and evaluation of the real-time buffet data are described, followed in a subsequent 
section by the spectral analysis of the same data. 

The manner by which the transonic maneuver tests were carried out is described 
below. Prior to entering the basic maneuver (the wind-up turn),- the pilot maintained 
the aircraft at the specified altitude and Mach number. The maneuver was executed by 
a combined turn and roll motion at maximum thrust. The aircraft was allowed to lose 
altitude so that the Mach number was maintained. In a typical wind-up turn at the test 
altitude range, 7,772 to 10,668 m (25,500 to 35,000 ft), the lost altitude was in the range 
of 152-305 m (500-1000 feet). During this time, the angle of attack was increased 
gradually to a maximum value. As the angle of attack increased, the buffet onset was 
first encountered. The buffet onset could be detected either based on instrumentation 
(e.g. , accelerometer under the pilot seat) or based on pilot perception. As the angle 
of attack reached its maximum, sustained structural vibration of the aircraft (i.e. , 
wing rock) took place. Wing rock usually was accompanied by severe rigid body 
oscillations such as the yaw and pitch motions. The structural vibrations and the 
rigid body oscillations both affected and degraded the aircraft's tracking ability. To 
terminate the maneuver, the pilot pushed the stick forward and returned the aircraft 
to a level position. 

Typical oscillographs giving time histories of the angle of attack, the sideslip 
angle, the roll angle, and the lift coefficient (computed) of a wind-up 
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turn are shown in Figures 11 and 12. The data were collected in a previous program 
(Reference 26) with M q = 0.89 and altitude h = 10,668 m (35,000 ft). The oscillographs 
clearly showed the oscillatory motions of the aircraft during the maneuver as the angle 
of attack increased. They are typical and similar to data collected in the present 
program and are described in the following subsections. 


3.1 PCM Data 

During the buffet flight test program, both the PCM and FM recording systems 
were activated. The PCM system recorded most flight condition parameters as well 
as the dynamic pressure data on the bottom surface of the right wing (4 transducers) 
plus Pressure Transducer No. 4 on the top surface (see Figure 3). The data passed 
through the presampling (low-pass) filters which retained the steady-state* and low- 
frequency components up to 30 Hz (120 Hz for super-commutated PCM data). The 
data were valuable in assessing the general behavior and response of the aircraft 
during the maneuver. This was in contrast to the FM data which were screened by the 
high-pass filters so that the dynamic components related to shock and shock oscillation 
could be amplified and studied in detail by eliminating the steady-state data. Typical 
PCM data recorded in Flight 825 are presented in this subsection. In Rim Number 5 
of the flight, the Mach number was 0.925 and the initial altitude was 10,668 m (35,000 ft). 
The leading edge and trailing edge flaps were in retracted position during this run. 

Figure 13 shows the pitch, roll, and yaw rate time histories during the wind-up 
turn. The elapsed time in seconds as recorded on the flight recorder magnetic tape 
is entered on top of the oscillographs. Also attached is the instantaneous angle of 
attack data retrieved from another PCM channel. The data showed that the aircraft 
entered the maneuver with a combined roll and yaw motion at least 1 second before the 
angle of attack increased significantly. At the time interval when the angle of attack 
was highest, significant low frequency roll oscillation was observed which was consis- 
tent with the general definition of "wing rock." 

Figure 14 gave the pilot seat lateral and normal accelerations Ny, N z during the 
maneuver. Referring to the figure, the peak load factor appeared at the time span 


*The steady-state data as defined here indicate the instantaneous mean values after 
the dynamic components above approximately 1/2 Hz have been extracted. Since the 
transonic maneuver was dynamic in nature, all the physical data varied with time, 
including the so-called steady-state data. 
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when wing rock took place. A negative load factor was recorded during the recovery 
stage of the wind-up turn. 

Typical PCM recordings of the dynamic pressure transducer "steady state" output 
are given in Figure 15 made for the same run (Run 5, Flight 825). The oscillographs 
were transcribed from a single PCM data set as compared to the digital tape data 
which were compiled from the super-commutated PCM data sets. (Reference 
Section 2.2.) Pressure Transducer Number 4 was located at 85% semi-span and 90% 
chordwise position at the top surface of the right wing. The data showed a substantial 
decrease in pressure (increase in lift) as the angle of attack increased. This is con- 
sistent with the flow behavior prior to and during the development of a shock induced 
separation wake on the airfoil surface. 

The remaining traces of Figure 15 give the steady-state and low-frequency 
pressure data on W. S. 128.31 (85% semi-span) at the bottom surface of the right wing. 

The three traces were for 20%, 60%, and 90% chordwise positions (Transducer Numbers 
19, 20, and 21) respectively. As the angle of attack increased, the pressure increased 
correspondingly at 20% and 60% positions, with the increase at the 20% position sub- 
stantially higher than that at 60%. A decrease in steady-state pressure was registered 
at 90% chordwise position due to the wake flow developed at the trailing edge area and 
the upstream shock on the top surface. The pressure data of Figure 15 were pressures 
relative to a static source registered in the manifold chamber at F.S. 349.5. At this 
altitude, 10,668 m (35,000 ft), the aircraft lost approximately 259 m (850 ft) during the 
maneuver. The reference pressure increased approximately 1372 N/M (0.199 psi) 
corresponding to the drop in altitude. This substantial variation in reference pressure 
during the transonic maneuver should be considered in evaluating the steady-state 
pressure data. 

Figure 16, taken from Reference 13, gives a qualitative description of the shock 
and separated flow development on a two-dimensional airfoil as the angle of attack is 
increased, while the Mach number is kept constant. Recognizing the significant 
difference of a two dimensional flow and the effect of the airfoil thickness, etc., 

Figure 16 is nevertheless included here to assist in visualizing the shock formation 
and movement and the corresponding steady state pressure changes. 

The set of PCM data as presented, based on Run 5, is concluded with Figure 17, 
which gives the normal and longitudinal accelerations (N^., N^) at the C.G. of the aircraft 
as well as the normal accelerations recorded at two stations of the right-hand wingtip. 
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The normal acceleration at the C.G. was similar in general pattern to those under the 
pilot seat (Figure 14, top trace), while the amplitude was somewhat lower at the eg, 
as was to be expected. The wingtip normal accelerations indicated that structural 
vibrations were superimposed over the rigid body motions of the aircraft. The dynamic 
components of the wing accelerations will be discussed later in the report. 

At lower altitude, 7,772 m (25,500 ft), and lower Mach number, M q = 0.75, a 
test run (Run 7) was made in Flight 825. In general, the steady-state pressure pattern 
was similar to the previous case (M q = 0.925, h = 10,668 m), while the amplitudes were 
somewhat lower for the lower Mach number case. Correspondingly, the load factor 
oscillographs showed the same general pattern with relatively lower peak amplitudes. 
The normal and longitudinal accelerations at eg and the wingtip normal acceleration 
data for Run 7 are given in Figure 18. The data may be compared with those of 
M q = 0.925 (Figure 17) to observe the Mach number effect. 









FLIGHT 657 FLIGHT TEST & BUFFET CRITERIA 

1 80.5 92.0 0 81 RUN 7 35KFT 0.89M WIND-UP TURN FLAPS 0/0 
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Figure 12. Lateral-Directional Oscillation Histories of a Transonic Wind-Up Turn, 
M q = 0.89, h = 10,668 m (35,000 feet) d = 0°, <5 f = 0° 


36 





0.925, h = 10,668 m 







Histories of Flight 825, Run 5, M - 0.925, h 10,668 m, 
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SHOCKWAVE: 
NO SEPARATION 



Figure 16. Successive Stages in an Incidence Increase at Fixed M Q (Pressure 
Distribution is Shown in Top Charts; Corresponding Flow Patterns 
Are Shown in Lower Sketches . . . From Reference 13) 
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Figure 17. CG and Right Wing Tip (Forward and Aft) Acceleration Time Histories 
of Flight 825, Run 5, M = 0. 925, h = 10, 668 m, 6 = 0°, 6 f = 0° 



3.2 FM Data 


In the buffet flight tests, 31 channels of data were recorded in the FM system 
in addition to the event marker. In order to obtain high-quality transient dynamic 
data, high -pass filters with a cutoff frequency of 1/2 Hz were used to eliminate the 
so-called steady-state components. The data list was given in Table 2 of the report. 
Specifically, of the 31 channels, 20 channels were used to record the dynamic 
pressure data acquired by the pressure transducers on the right wing top surface 
(Figure 3). Six channels were assigned to strain gage inputs. The strain gage data 
were used to compute the wing structural loads. One channel was used to record 
the right hand aileron hinge moment, and the remaining four were for dynamic loads 
data recorded by accelerometers. 

A set of typical FM data acquired in Flight 871, Run 2 (M q = 0.925, h = 10,668m 
or 35,000 ft) is presented in this subsection. During Run 2, the leading edge and 
trailing edge flaps were extended to<5 n /«i = (4°/l2°). In general, the dynamic pressure 
data at the outboard wing stations were less significantly affected by the extension of 
the flaps as compared to the lower Mach number case. The corresponding pressure 
data at the inboard stations, especially those wing stations covered by the trailing 
edge flaps (see Figures 1 and 3), were to a larger degree influenced by the flap 
settings. 

The dynamic pressure data acquired by the FM system are presented in Figures 
19 through 25. Each figure has three oscillographs that are properly identified 
(see Figure 3 for station number identifications). On top of these figures, the time 
counts in seconds are noted, as well as the instantaneous angle of attack data. The 
traces covered 14 seconds during which a wind-up turn was performed. As the angle 
of attack was increased, the first major noticeable pressure oscillation appeared at 
t = 073. 0 at station Number 2 (Figure 19). The pattern indicated oscillation of the 
shock front as a result of downstream separation during the maneuver. Downstream 
separation can be seen in the pressure fluctuation at stations 3 and 4 at t = 073. 0 This 
can be defined as the buffet onset point. It seemed that the pressure drop across the 
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shock front may be determined, at least approximately, by the amplitude of the pres- 
sure oscillation (about 6000 N/M or 0.87 psi). The pressure oscillations defining 
the location of the shock induced separation did not show significant influence on the 
fluctuating pressure at upstream locations. They did affect the downstream stations 
where the separated flow took place (see Pressure Traces Numbers 3, 4, 22). As the 
angle of attack was further increased (to approximately 6.2° at t = 076. 0, the separa- 
tion boundary (as defined by the shock front oscillation) moved to Station No. 1 (20% 
chordwise position) at 85% semi-span section. A comparison of Pressure Traces 
Numbers 1 and 2 indicated that as the angle of attack was increased, the shock-induced 
separation boundary stayed in the neighborhood of Station Number 2 for approximately 
two seconds, during which a increased 2.2°. It moved more speedily through Station 
Number 1 toward the leading edge of the wing. The pressure traces behaved in a some- 
what random fashion when the angle of attack continued to increase and the local station 
was under the separated flow. 

Similar shock and pressure fluctuation behavior was observed at W.S. 110. 30 
(72.8% semi-span, Pressure Station Numbers 5, 6, 7, 8, 9) and W.S. 93.0 (61% semi- 
span, Pressure Station Numbers 10, 11, 12). A special point of interest was noted in 
Pressure Trace Numbers 6 and 11 where the distinct pressure oscillation appeared 
and disappeared in less than 0.20 second. The pressure across the shock, as indicated 

by Pressure Trace Numbers 5 and 10, was almost identical to those based on Pressure 

2 

Trace Number 1 (6000 N/M ). The approximate time when the complete top wing 
surface was under separated flow was identified through Pressure Trace Number 15 
when large amplitude random oscillation appeared abruptly at t = 078.9, with a =9.4°. 
Pressure Stations 14 and 18 were on the trailing edge flap surface. Except for a 
somewhat higher peak-to-peak oscillation amplitude, no major variations of these 
pressure traces were observed as compared to those of nearby stations. 

Typical wing structure strain gage data were recorded as the third trace of 
Figure 25, three traces of Figure 26, and the first and second traces of Figure 27. 
These strain gage data were used to determine the dynamic loads at a given section 
based on empirical formulas established in ground tests. Figure 28 gives the right- 
hand aileron hinge moment and normal accelerations at eg and the right-hand wingtip 
(forward position). The aileron hinge moment showed fairly high-frequency oscillations 
which may be damped out by the control linkage. Also, good correlation was observed 
between the aileron hinge moment data and the pressure data picked up at a location 
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close to the hinge (Pressure Number 12). The accelerometer data for N z at the C.G. 
and the wingtip during buffet onset and wing rock will be discussed later in the section 
dealing with analytical correlation. 
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Figure 19. Oscillographs of Pressure Station Numbers 1, 2, 3, Recorded in Run 2, 
Flight 871. M = 0.925, h = 10,668 m, 6 =4°,«, = 12°. 
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Flight 871. M = 0.925, h = 10,668 m 
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Figure 21. Oscillographs of Pressure Station Numbers 7, 8, 9, Recorded in 
Run 2, Flight 871. M = 0. 925 , h = 10,668 m , 6 =4°,i f =12°. 
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Figure 24. Oscillographs of Pressure Station Numbers 16, 17, 18, Recorded 
in Run 2, Flight 871. M = 0. 925 , h = 10, 668 m , a = 4° , <5 f = 12 
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Figure 25. Oscillographs of Pressure Station Numbers 22, 23, and Gage RW32B18 
Recorded in Run 2, Flight 871. M = 0.925, h = 10,668 m, 6 = 4°, 
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in Run 2, Flight 871. M = 0.925, h = 10,668 m, 
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Figure 27. Oscillographs of Gages RW68B18, RW68T30, Recorded in Run 2, 
Flight 871. M = 0.925, h = 10,668 m, 6 = 4 °, 6 . = 12 °. 
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3. 3 Flap Setting and Mach Number Effects 

Through visual examination of the real-time pressure data corresponding to 
identical Mach number and altitude , but with varying leading edge and trailing edge 
flap settings, the effect of the latter to the shock pressure oscillation patterns may be 
evaluated. Other flap setting effects such as the development of the separated flow 
region with increasing angle of attack, the influence on dynamic loads, etc. , will be 
discussed later in the report. 

Figures 29, 30 give the pressure transducer oscillographs obtained at Station 
Numbers 1, 2, 3, and Numbers 10, 11, 12 respectively. The oscillographs were 
obtained based on Flight 825, Run 5, where M q =0.925, h = 10,668m (35,000 ft) and 
both leading edge and trailing edge flaps were retracted. The data may be compared 
with Figures 19, 22 respectively, where the Mach number and altitude were identical 
while the flaps were set at (4°/l2°). It may be noted that the flap retracted maneuver 
was conducted at a faster pace. The elapsed time of the maneuver lasted approximately 
seven seconds so that the peak angle of attack (14. 4°) was reached at a comparatively 
short time after the first shock appeared on the wing surface. 
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Flight 825. M = 0. 925, h = 10,668 m 
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Figure 30. Oscillographs of Pressure Station Numbers 10, 11, 12, Recorded in 
Run 5, Flight 825. M = 0. 925, h = 10,668 m, <5 =0°,<5 f =0°. 



The Mach number effect may be observed in a similar manner by examining the 
real-time data obtained in two runs where other conditions were equal except for 
Mach numbers and flight altitudes with the q value unchanged. 

Figures 31, 32 are the fluctuating pressure tracks at Stations Numbers 1, 2, 3, 
and Numbers 10, 11, 12 respectively for Run 12 of Flight number 871. The run was 
executed at M q =0.75 and h = 7, 772 m (25, 500 ft). The leading edge and trailing edge 
flap settings were (4°/0°). At time count 538, a brief reduction of the angle of attack 
was noted. In reviewing the other channels of the flight test data, it was determined 
that the temporary reduction in the angle of attack was caused by the pilot. A com- 
parison of these figures with Figures 19, 22, (M q =0.925, h = 10,668 m) revealed the 
difference in buffet pressure behavior at a lower Mach number. The difference in 
trailing edge deflection is not considered to have significantly affected the comparison 

Figure 33 gives the pressure traces at Stations 1, 2, 3, for M q = 0.75, 
h = 7,772 m while the leading edge and trailing edge flaps were retracted. These 
data may be compared to those of Figure 31 where the leading edge flap was extended 
(<5 n /c5f = 4°/0°). Figure 31 showed that the pressure oscillations at Station Numbers 
1 , 2 were somewhat delayed as the angle of attack was increasing when compared to 
Figure 33. This delay in buffet onset was attributed to the deflected leading edge 
flap (4°). 
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Figure 31. Oscillographs of Pressure Station Numbers 1, 2, 3, Recorded in 
Run 12, Flight 871. M = .75, h = 7,772 m, 6 =4°,<5 f =0°. 
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Run 12, Flight 871. M = .75, h = 7,772 m, 
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Figure 33. Oscillographs of Pressure Station Numbers 1,2, and 3, Recorded in 
Run 7 , Flight 825 , M = .75, h = 7,772 m,d =0°,<5 f = 0 o . 



3.4 Additional Observations 


In this subsection, the real-time data acquired from three runs in Flight 825 
are further reviewed. Flight 825 was conducted with leading edge and trailing edge 
flaps retracted (3^/ij ~ 0°/0°). The runs (Numbers 7, 6, 5 according to the sequence 
of presentation in this subsection) were conducted at M q =0.75, 0.85, 0.925 respectively. 

Figure 34 illustrates schematically the variation of the separation boundary with 
time obtained from the real-time pressure data such as presented in Figures 29, 30, 

33. In determining the boundary details where the local fluctuating pressure data 
were unavailable, use was made of the tuft data obtained in a previous flight test pro- 
gram (Reference 26). Figure 35, 36 presents the variation of peak-to-peak pressure 
turbulence with time and a for transducers located along three spanwise station (W.S. 
128.31, 93.0, 50.00). The values, measured from time histories, were for M q = 

0.75 and 0.925 (Run Numbers 7, 5) respectively. 
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(a) RUN 7, M 0 = 0.75, h = 7,772 m 



(b) RUN 6, M 0 = 0.85, h = 9,449 m 



(c) RUN 5, M 0 = 0.925, h = 10,668 m 

Figure 34. Upper Wing Surface Separation Boundary Mapping for Three Runs of 

Flight 825, (a n /a f = 0°/0°) 
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TRANSDUCER NO. 



Figure 35. Peak-To-Peak Pressure Histories at Various Wing Stations, 
Flight 825, M = 0. 75, h = 7, 772 m, 6 = 0°, <5 = 0° 












Referring to Figure 34(a) corresponding to a maneuver conducted at M q =0.75, 
the buffet onset appeared to be caused by leading-edge-induced separation. The regions 
of separation developed at the outboard wing first and spread inboard. From the initial 
separation at the wingtip to the complete separation over the entire wing, a time lag of 
approximately 2. 2 seconds was detected. During this period the aircraft a increased 
from approximately 3.7° to 11.5°. Since the separation was leading-edge-induced, 
there was not much chord-wise variation in the dynamic pressures with time at a fixed 
span-wise station. Once separation occurred at the leading edge, the entire region 
downstream was affected. 

The pressure turbulence along a given chord appeared to peak shortly after the 
onset point with a subsequent decrease for higher angle of attack (see Figure 35). This 
was to be expected for leading-edge-induced separation; i.e., at the critical a or slightly 
greater, the separated shear layer was relatively close to the surface. Therefore, the 
resulting turbulence was impressed on the wing surface. With increasing a, the shear 
layer moved away from the surface. The resulting turbulence was attenuated somewhat 
before it was felt at the surface. 

The separation boundary movement and peak-to-peak pressure histories for 
M q = 0.925 were given in Figures 34(c) and 36. For this case, the buffet onset appeared 
to be caused by shock-induced separation. The growth of the separation region was 
characterized by a movement in the chordwise (forward) direction especially at the 
early stage of the maneuver. The shock also spread from the wingtip to the inboard 
regions of the wing as the angle of attack increased. 

Referring to Figure 36, the high values at the buffet onset point at W.S. 128. 3 
and 93. 0 could be attributed to the shock oscillating over the particular transducer. 

The inboard pressure (W.S. 50.0) did not show any strong shock ocillations as at the 
other two wing stations. The separation in this region appeared to develop simulta- 
neously similar to the M q = 0. 75 data. In this case, the inboard separation resulted 
from the inboard spread of the outboard separation region. It is interesting to note 
that once the flow on the top surface of the wing was completely separated, the level 
of the peak-to-peak turbulence was similar to the leading edge induced data at 
M q = 0. 75. This was to be expected for as the shock approached the leading edge, the 
resulting downstream separation could be expected to resemble leading-edge-induced 
separation characteristics. 

The development of the separation boundary for M q =0.85 (Figure 34 (b)) followed 
a pattern which was intermediate in nature and resembled to some extent either of the 
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two extreme cases (M Q =0.75, 0.925). The separated flow seemed to be induced ini- 
tially by a weak shock (i.e. , a weaker shock as compared to the M q = 0. 925 case). 

As the angle of attack increased, the shock moved toward the leading edge swiftly. 

In the later phase of the maneuver, the expansion, of the separation boundary in the 
inboard portion of the wing was very similar to those of the M q = 0. 925 case. 

Figure 37, taken from Reference 26, indicates the separation boundaries of the 
F-5A Wing based on tuft data. The data were collected at h = 10, 668 m (35, 000 ft) 
corresponding to three Mach numbers (M q = 0.81, 0.86, 0.91). During the flight, 
the leading edge and trailing edge flaps were retracted. A major difference in the 
test aircraft configuration on which the data of Figure 37 were based was that 
the aircraft had a clean wing, while wingtip missiles were attached in the present 
program. Referring to Figure 37, the case corresponding to M q =0.81 indicated 
that the separation was leading-edge-induced, as was the case for Run 7, Flight 825 
(M q = 0. 75, h = 7, 772m) discussed previously. For the remaining two cases of 
Figure 37 (M q =0.86, 0.91), the separation seemed to be induced by local shock. 


INITIAL FLOW SEPARATION 
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Figure 37. Wing Tuft Results, h = 10,668 m, <5 = 0°, <5 f = 0° (From Reference 26) 











Section 4 


SPECTRAL AND STATISTICAL PROCESSINGS 
OF BUFFET DATA . 


In the previous section, the real-time buffet pressure and response data were 
presented. The dynamic pressure data covered a wide frequency range. Excluding 
distinct physical behaviors such as the shock pressure oscillations, etc. , the data 
within a given time span may be considered random. Assuming that the random 
data were stationary within a limited time span, standard spectral processing tech- 
niques may be applied. In this manner, certain general conclusions may be drawn 
based on the processed pressure and response data. Furthermore, the spectral 
processing made analytical correlation feasible in case the basic aircraft response 
function was assumed linear. 

The spectral processing consisted of the determination of the power spectral 
densities of the pressure and response data, as well as the cross-spectral densities 
of any two sets of data. For the latter case, the spatial correlation of the data sets 
came into play. In general, a high degree of spatial correlation may be observed if 
the two sets of data are caused by a single source, or if a given physical phenomenon 
is moving in space at a distinct speed (convection speed). In addition to spectral 
processing, other statistical processing techniques were also applied to the buffet 
data acquired in the flight test program . 


4.1 Processing Techniques 

The spectral and statistical processings of the digitized buffet data applied in 
the present program followed the common procedures used in most flight test and 
other experimental programs. For completeness in presentation, and to prevent 
ambiguity in definitions of certain functions, a brief outline of the processing tech- 
niques is described below, to be followed by systematic presentation of the processed 
results. 
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Consider a stationary ergodic random process x(t). Let x(t) be so truncated 
that it becomes zero outside the interval (o, T). The truncated function is denoted 


by x T (t). The Fourier transform of x T (t) exists provided that the absolute value of 
x T (t) is integrable and x T (t) is bounded: 

x (t) = /* F (f) e 12rft df 

T' -oo X 

(4-1) 

where F (f) -"f / T x (t) e 

x ' J o T 

(4-2) 

jjj afc 

If F (f) denotes the complex conjugate, then F (-f) = F 

X XX 

(f), since x T (t) is real in the 

present case. 

The auto-correlation function is defined as: 

- 

R (r) = lim y/ T x(t) x(t + f ) dt 
x T-*°° 1 ° 

(4-3) 

and R x (-r) = R x (r) , R x (o) > | R x (t) | for all t 

(4-4) 

For stationary data, the power spectral density function and the auto-correlation 
function are related by a Fourier transform as follows: 

0 (f)= 2^11 l(T )e" i2rfT dr 

A -CO A 

(4_5) 

= 4/ C0 R ( T ) cos 2irfrdT 

® X 

The mean square value of x(t) is given by: 

■ — ' T 

x 2 = lim j f x 2 (t) dt 

T — *oo 0 

(4-6) 

■ C '’x « df = R x (°) 


Applying Parseval's theorem, the power spectral density function may be expressed 
as follows (see Reference 47): 

0 x (f ) = 2T | F x( f )| 

(4-7) 

Similarly, for two sets of stationary ergodic random data 

x(t) and y(t), the cross- 

correlation function is defined as: 

. 

R v (T) = lim J f T X (t) y (t +t) dt 
T— *co o 

(4-8) 

and R r (-A = R (t) 

xy ' ’ yx ' ’ 

(4-9) 
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The cross -spectral-density function and the cross correlation function are related by 
the Fourier transform: 


^ (f) = 2 _C R xy (T>e ' i! ” fTdT 

Two useful formulations relating auto- and cross-power spectra are: 


(4-10) 


(4-11) 


0 X y ( f ) 2 <0 x (f)0 y (f) 


(4-12) 


The coherence function is a real -valued quantity given by: 

y 2 x 1 

xy w 0 x (f) 0 y (f) 1 


(4-13) 


When y (f) - 0 at a particular frequency, x (t) and y (t) are said to be uncorrelated 

at that frequency. If x (t) and y (t) are statistically independent, then Y^ (f) = 0 for 

all frequencies. When y 2 (f) = 1 for all f , then x (t) and y (t) are said to be fully 

xy 


coherent or perfectly correlated. A quantity equivalent to coherence function in the 
time -domain is the correlation function coefficient (normalized cross-covariance 
function) which is defined by: 


R xy (T) - "y 

P (T)= ===== = = T . ■ 

-1 < P^ (r) < 1 for all r's 


(4-14) 

(4-15) 


where and p are the mean values of x (t) and y (t) respectively. The function p (r) 
y 

measures the degree of linear dependence between x (t) and y (t) for a displacement 
of r in y (t) relative to x (t) . 

The cross-correlation measurement between two sets of random data may yield the 
time delay information when the data are properly correlated. For this case, the 
time delay may be determined by noting the time displacement associated with an 
observed peak in the cross-correlogram. 

The probability density function for random data x (t) describes the probability that the 
data will assume a value within some defined range x< x(t) < x + A x at any time instant: 
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(4-16) 


p (x) = lim lim 1 T x 
Ax— o T— co "T 

where T is the total amount of time that x(t) falls inside the range (x,x+ Ax) during 
an observation time T. The cumulative probability distribution function which 
describes the probability that the instantaneous value x (t) is less than or equal to 
some value x is defined by 

P(x) = Prob [x(t) < x] = / X p(|) d£ (4-17) 

- oo 

and it is bounded by zero and unity. 

4.2 Digitized Power Spectra 


The digitizing process converted the continuous flight test data into discrete 
numerical data. The sampling rate for high-frequency digitization procedure used 
in the subject program was 5,000 samples per second, corresponding to a cut-off 
frequency of 2,500 Hz. The highest frequency in the presented data was limited to 
1,000 Hz, which was considered sufficient to handle the buffet pressure and response 
data of a full scale aircraft. In order to achieve detail spectral coverage in the low- 
frequency range where all major structural modes of the aircraft were located, a 
separate low-frequency processing procedure was used. 


The sampling rate for low-frequency digitization procedure was 1, 000 samples 
per second. The frame rate was 0.001 second and the cutoff frequency was 500 Hz. 
Sampling rate other than the ones quoted above were also used depending on specific 
application. In order to avoid the aliasing problem, the high-frequency components 
were filtered out prior to sampling. 


There are two methods to generate the power spectral data: 

1. Fourier transform method — The digital time history functions are first 
converted to Fourier transforms utilizing Equations (4-1) and (4-2) as 
follows: 
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where N is the sample size, and At is the time increment (equal to a 
multiple of the frame rate) which has been selected to produce a sufficiently 
high-frequency cutoff, f,, = 1/(2 At). The associated basic time span T is 
N At. The frequency increment is l/T. 

Application of equation (4-7) yields the power spectral data. For high 
frequency data, more time segments have to be used. The final results 
are obtained by averaging over all the segments. 

2. Correlation method — The digital data are first converted to auto correla- 
tion functions utilizing Equation (4-3) as follows: 


R / = R (iA t) 




n = 1 


x n x n +i 


f= 0, 1, 2,...,L (4-20) 


where L is the maximum correlation lag. Let B g be the desired equivalent 
resolution bandwidth for PSD calculation, then 

L - -g-ij- • (4-21) 
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The sample size N is chosen such that (Reference 47) 
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where « is the normalized standard error for spectral calculations, and the associated 
minimum time span T is T = N 4t. Equation (4-22) establishes the normalized stan- 
dard error of power spectra as a function of the resolution band width B g and the 
sample data time span T. This equation places certain restrictions on the spectral 
processing of the transonic maneuver data where the sampling time is limited. This 
point will be discussed further in a later subsection. 


In digital processing of the power spectral data, Reference 47 suggests that 
the raw estimate 0(f) of a true power spectral density function 0 X (f) be calculated 
only at L + 1 discrete frequencies: 

kf , 

f= -L 5 ‘ = 2LAt k = 0,1,2 L (4-23) 

At these discrete frequency points, 
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R^ cos (-^ ) + (-l) k R l ] 


(4-24) 




A final estimate of the power spectral density 0^ is obtained using the following 
smoothing process. Smoothing is necessary since the raw estimate given by Equation 
(4-24) is an inefficient estimate of the true spectral density (Reference 47). 


0 O = 0. 5 0 Q .+ 0. 5 
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(4-25) 
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In digital computation, at the discrete frequencies f = k/(2L4t), k = 0, 1, 2. 
L , the coherence function is estimated by: •„ 


7 (k) = 

xy ' ' 


0xy (k) 


0 X (k) 0 y (k) 


(4-26) 
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4.3 Typical Pressure and Response Power Spectra 


In this subsection typical power spectra obtained from the real-time buffet data 
are presented. The variation in spectral make-up trends corresponding to various 
phases of the maneuver is discussed. (A systematic presentation of the bulk of the 
spectral data is given in subsection 4.4.) Also included in this subsection are the low- 
frequency spectral data which are important in aircraft response studies. The low- 
frequency response is significant because all the major structural modes are located 
in this frequency range. 

The pressure data were processed and converted into power spectra for a wind- 
up turn at M q = 0. 75 and h = 7, 772 m, Run 7, Flight No. 825. Figure 38 shows the 
pressure power spectra obtained at the top surface of Wing Station 128.31, 85% semi- 
span and 90% chord position for a separated flow condition. Pressure Station No. 4, 
Figure 3. The power spectrum in the frequency range above 100 Hz seems to follow 
the -5/3 slope indicated by the inserted broken line. While not claiming any physical 
similiarity between the separated flow behavior described here and the case of the 
homogeneous turbulent flow, it is worthy of noting that the -5/3 slope is predicted by 
the theory of universal equilibrium in homogeneous turbulent flow. (See Reference 38, 
"The Structure of Turbulent Shear Flow," by A. A. Townsend, page 45.) The time span 
during which the spectral data of Figure 38 were acquired was t = 363.4 - 364.22. 
During this time, the flow on the complete upper wing surface became separated, as 
was shown in Figure 34(a) . A typical power spectrum plot for an inboard pressure 
station located on the trailing edge flap, Pressure Station No. 18, Figure 3, acquired 
at the same time interval is shown in Figure 39. Figure 39 shows that the PSD level 
was quite high at this time; it represents an overall peak in spectral power for Station 
No. 18 during the complete maneuver. The high spectral power at Station No. 18 was 
attributed to the high turbulence when the local flow became separated where the shear 
layer was relatively close to the surface. 

In order to examine the variations of the wing surface pressure power spectrum 
with angle of attack, the mean PSD function for various angle of attack ranges encoun- 
tered in the wind-up turn maneuver was calculated. (The term "mean PSD" was used 
because a finite time increment corresponding to the angle of attack range variation 
was used to define the PSD function; hence, the requirement of stationarity was not 
exactly satisfied.) Furthermore, very narrowband spectral peaks were dropped to 
reveal the overall pattern. 
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Typical results of the type described above for Pressure Station No. 4 are shown 

in Figure 40. In Figure 40 the PSD function for M q = 0. 75 is presented for five angle 

of attack ranges. For a = (2.2° - 3.8°), the flow was essentially unseparated on the top 

2 2 

of the wing. Except for very low frequencies, a PSD in the 47.54 (N/M ) /Hz 
(10 psi /Hz) range was observed. The pressure PSD showed a substantial increase 
as the angle of attack was increased to the range (4. 3° - 6. 0°). For the remaining 
period of the maneuver, the PSD below 100 Hz did not vary greatly. Above 100 Hz, 
some amplitude changes were observed. The drop in PSD with an increasing a in this 
frequency range is believed due to the thickening of the separated flow region above the 
station where the pressure was measured, as mentioned previously in the report. 

Two similar PSD trend plots for Pressure Station Nos. 4 and 1 are given in 
Figures 41 and 42 corresponding to M q = 0.925, h = 10,668 m. Figure 41 shows the 
relatively small change in spectral power level at Station No. 4 throughout the tran- 
sonic maneuver across the frequency band. Compared to Figure 40, it shows clearly 
the Mach number effect at the early stage of the maneuver. For M q = 0.75, the flow 
at W.S. 128.31 did not become separated until the angle of attack z’eached at least 4° 
or more. Very low spectral level was thus recorded in Figure 40 corresponding to 
a= (2.2° - 3.8°). For M q = 0.925, the flow became separated at the same section at 
a smaller angle of attack. This trend was clearly reflected in Figure 41 in the rela- 
tively high spectral level at a = (2.2° - 3.3°). Another point of interest was the simi- 
larity in spectral patterns for frequencies above 100 Hz corresponding with the two 
Mach numbers. 

When compared to Figure 41, the spectral plot for Pressure Station No. 1 (Fig- 
ure 42) shows the essential difference in dynamic pressure behavior at various chord- 
wise positions of a wing section at M q = 0.925. At low angle of attack, Station No. 1 
registered very-low-level spectral power because the local flow was unseparated and 
was unaffected by the shock which appeared downstream. However, part of the low 
level spectra may be attributed to the instrumentation noise. 

The PSD of Pressure Station No. 1 (Figure 42) was highest at a = (4.2° - 6. 2°) 
when the local shock caused an oscillating pressure pattern. As the maneuver pro- 
gressed, the pressure PSD level dropped somewhat. 
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In Figures 40-42, the PSDs corresponding to the recovery phase of the maneu- 
vers after the peak angle of attack has been reached are not included. In general, 
the PSD levels corresponding to the recovery phase were comparatively high, which 
were attributed to the fluctuating pressures generated by the wing during the rapid 
recovery as well as the local fluctuating pressures generated by the deflected control 
surfaces. Sample data corresponding to this phase are given in the following sub- 
section. 
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Figure 38. Power Spectrum of Buffet Pressure at Wing Station 128.31, 90% Chord 
Position, for M = 0. 75, h = 7, 772 m, 6 =0°, 6 C = 0° 
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Figure 39. Power Spectrum of Buffet Pressure at Wing Station 50. 00, 80% Chord 
Position, for M q = 0. 75, h = 7,772 m, <5 n = 0°, <5^=0° 
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Figure 40. Trends of Buffet Pressure Spectra at Station No. 4 Corresponding to 

Various Angles of Attack During a Wind-Up Turn, M = 0.75, h = 7, 772 m 
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4.4 Power Spectra Corresponding to Various Phases of the Transonic Maneuver 


In this subsection, auto power spectra of pressure and responses are presented 
in a systematic manner corresponding to various phases of the transonic maneuver. 

As mentioned previously, the transonic maneuver was dynamic in nature. The local 
buffet pressure and the response data were nonstationary in a strict sense. Within a 
short time interval (1 or 2 seconds) when no drastic change in local physical behavior 
took place, the data may be considered random and near-stationary so that the spectral 
processing technique may be applied. 

For data acquired in Run 5, Flight 825, (M q = 0. 925, h = 10, 668 m, 4^/^ = 0°/0°), 
five time segments were chosen for spectral processing. Each time segment repre- 
sented 1.025 seconds. Roughly, the five time segments may be classified as follows 
(also see Figure 34(c)). 


Designation 

Starting Time 

a Initial 

Description 

A 

334.0 

m 

Initiated wind-up turn. Shock appeared 



■ 

at localized area. 

B 

335.03 


Buffet onset. 

C 

336.06 

7.8° 

Separation region expanded. 

D 

. 

337.09 

12.3° 

Separated flow covered the complete 
wing surface. 

E 

338.12 

13.9° 

Recovery initiated. 


The five PSD plots for each function corresponding to the time segments (A) 
through (E) are presented in one figure. For Run 5, seven pressure data sets were 
selected. These are Pressure Station Numbers 1, 3, 4, 6, 12, 18, and 23, whose 
PSD's are shown in Figures 43-49 respectively. Following the pressure PSD plots are 
four sets of response PSD plots. These are normal acceleration PSD at C.G. (Figure 50), 
normal acceleration PSD at forward and aft stations of right-hand wingtip (Figures 51, 
52), and right-hand aileron hinge moment PSD (Figure 53). 

Of the PSD plots described above, it may be mentioned that starting at t = 336. 06, 
some low-frequency vibrations appeared in the aircraft. These vibrations were 
noted in the acceleration PSD plots as minor peaks. All the acceleration PSD plots 
(Figures 50-52) showed very high peaks between 400-500 Hz. It is possible that 
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these peaks were caused by the aircraft electrical noise because of the insufficient 
shielding of the lead wires and other unidentified disturbances. (The accelerometer 
harness was installed inside the aircraft perviously to acquire low- frequency 
response data.) 
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Figure 43. PSD's of Pressure No. 1 Based on Various Segments of 



Figure 44. PSD's of Pressure No. 3 Based on Various Segments of 
Real Time Data, M = 0. 925, h = 10, 668m, d/d f = 0°/0° 



Pressure No. 4 Based on Various Segments of 






Figure 46. PSD's of Pressure No. 6 Based on Various Segments of 
Real Time Data, M = 0. 925, h = 10, 668m, <5 /<5 f = 0°/0° 
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Figure 47. PSD's of Pressure No. 12 Based on Various Segments of 
Real Time Data, M = 0.925, h = 10, 668m, 6 / = 0°/0° 




Figure 48. PSD's of Pressure No. 23 Based on Various Segments of 
Real Time Data, M = 0. 925, h = 10, 668m, 6 /i. = 0°/0° 
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Figure 50. PSD's of C.G. Normal Acceleration Based on Various Segments of 
Real Time Data, M = 0.925, h = 10,668 m, 6 /6 f = 0°/0° 




Figure 51. PSD's of Righthand Wing Tip Acceleration Forward Based on Various Segments of 
Real Time Data, M = 0.925, h = 10,668 m, 6 /d r = 0°/0° 




Figure 52. PSD's of Righthand Wing Tip Acceleration Aft Based on Various Segments of 
Real Time Data, M = 0.925, h = 10,668 m, 6 /6 f = 0° /0° 
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Figure 53. PSD's of Right Hand Aileron Hinge Moment Based on Various 
Segments of Real Time Data, M = 0.925, h = 10, 668m, 6 /i f = 0°/0° 



Corresponding to the eleven sets of PSD plots for Run 5, Flight 825 (M q = 0.925) 
presented in Figures 43-53, five sets of PSD plots were generated based on Run 7 of 
the same flight (M q = 0. 75). The time segments for the latter run were chosen in such 
a manner that the angle of attack at the beginning of each time segment was comparable 
to those of Run 5, with the exception of the fifth time segment where the initial angle of 
attack was substantially lower for Run 7. The five sets of PSD data are shown in Fig- 
ures 54-58. Figures 54, 55, and 56 are for Pressure Stations No. 1, 3, and 23. Fig- 
ures 57, 58 show the PSD's for the normal eg acceleration and the right-hand wingtip 
acceleration, (forward station) respectively. 


A comparison of Figures 43 and 54 showed that for Pressure Station No. 1 (20% 
chord and 85% semi-span) the pressure PSD was substantially higher for M q = 0. 75 at 
the initial stage of the maneuver. This was due to the leading-edge-induced separation 
corresponding to M q = 0. 75 discussed previously in reference to Figure 42. For an 
inboard station such as Pressure Station No. 23, a comparison of Figures 48 and 56 
shows close similarity in the two sets of the PSD plots except in the low-frequency 
region. 

A typical set of fluctuating pressure PSD plots for Run 6 of Flight 825 at M q = 0. 85 
is shown in Figure 59. The PSD's are for Pressure Station No. 1. These plots may 
be evaluated against similar data for M q = 0. 925, 0. 75 presented above (Figures 43, 

54). In the case of Run 6, six (6) time segments were chosen covering the maneuver. 
The successive segments were designated A through F. Referring to Figure 59, the 
spectral trend plots reflected the relatively low fluctuating pressure prior to local 
flow separation (Segment A), the increase in power when the shock induced separation 
passed the locality (Segment C). The PSD level dropped with further a increase when 
the separation boundary moved inboard (Segment E). The behavior was similar to the 
case of M q = 0.925 (see Figure 42). Subsequently, the PSD level increased at the end 
of the maneuver (Segment F) 
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Figure 54. PDS's of Pressure No. 1, Based on Various Time Segments of 
the Real Time Data, M =0.75, h=7,772m, 6 /i f = 0°/0° 







Figure 55. PSD's of Pressure No. 3, Ba sed on Various Time Segments of 
the Real Time Data, M = 0. 75, h = 7, 772m, 6 /d f = 0°/0° 
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Figure 58. PSD's of Righthand Wing Tip Acceleration (Forward) Based on Various Time Segments of 
Real Time Data, M = 0.75, h = 7,772 m, 6 / 6 . = 0°/0° 







Figure 59. PSD's of Pressure No. 1 Based on Various Segments of Real Time Data. 

M = 0.85, h = 9,449 m, <5 /tf f = 0°/0° 





4.5 Correlation and Coherence Study of Buffet Data 

Following the approach in Subsections 4.1 and 4.2, the cross-correlation 

functions for two sets of digital data x^ and at lag numbers l = 0, 1,2 ..., L may 

be defined as R (l At) based on Equation (4-8). Furthermore, the cross-correlation 

— 

coefficient p and the raw estimate of the cross spectral density 0^ may be defined 
based on Equations (4-14) and (4-10) respectively. A smoothed estimate for the cross- 
spectral density 0 is obtainable using the same smoothing process of Equation 
(4-25). The coherence function is determined by Equation (4-26). 

Two types of cross-correlation and coherence functions are presented in this 
subsection: the first type is the correlations between wing root structural loads 
and the dynamic pressures; the second is the correlations between the pressure 
transducer pairs. All the data are taken from Flight 871, Run 2, low-frequency 
digital tape. The time span covers the period from 079. 5 to 081.6 second. The time 
increment is 0.002 second and the frequency increment is 0.488 Hz. The maximum 
time lag is 1.024 seconds. 

Three sets of data are presented covering the correlations between the wing root 
structural loads and the computed pressure loads. The computed pressure loads are: 

Total Bending Moment 
LPY = / py d A 

y = bending moment arm from W.S. 32 

Total Torsional Moment 
XPX =f pxd A 

x = torsional moment arm referred to elastic axis at W.S. 32 

Summation of the Measured Pressures 
£P = f pdA 

Corresponding to the computed pressure loads, correlation data were computed 
as indicated below. 
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1. The lower part of Figure 60 shows the correlation (as defined in Equation 
4-8) between the bending moment at W.S. 32. 0 and the total bending moment 
contributed by the measured buffet pressures. The maximum cross- 
correlation coefficient | p \ = 0. 573 occurs at time lag t = -0. 022 second. 

(See Equation 4-14 for definition of p.) 

2. Figure 61 shows the correlation between the torsional moment at W.S. 32.0 
and the total torsional moment contributed by the measured buffet pressures. 
The maximum cross-correlation coefficient j p ] = 0.4 occurs at r = -0.016 
second. 

3. Figure 62 shows the correlation between the shear force at W.S. 32.0 and the 
summation of the measured buffet pressures. The maximum cross-correla- 
tion coefficient | p| - 0.666 occurs at r = -0.014 second. 

For the above three cases, the maximum correlation coefficients appeared at a 
small negative r. This indicates that corresponding to the major modal responses, the 
buffet pressure loads always lead the structural loads by a small phase angle. In other 
words, the impedance of the aircraft (including structural and aerodynamic damping) is 
such that a small positive damping exists corresponding to the major response modes. 

The coherence function plots show large-amplitude changes from frequency to 
frequency. A contributing factor to the large amplitude fluctuations is the fact that the 
combination of a relatively small equivalent resolution frequency and the limited sample 
time span caused a large normalized mean square error of the spectral data. It is 
worthy of note that the coherence function between the pressure pairs showed less 
fluctuation when the same procedure was applied to the same batch of low frequency 
data. This point will be discussed in the following paragraphs. 
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Figure 60. Correlation Between the Wing Root Bending Moment and the Corresponding 
Moment of Pressures, M q =0.925, h = 10,668 m, <5 n =4°, =12°. 
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Figure 61. Correlation Between the Wing Root Torsion and The Corresponding 

Moment of Pressures, M =0.925, h = 10,668 m, 6 =4°, 6 r = 12°. 
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In the following, the correlation of the pressure data for two transducer locations 
is studied. Specifically, the cross correlation, the coherence function, as well as the 
phase angle of the' cross spectra of the pressure data yield information about the pres- 
sure spatial correlation. For instance, if a condition of pure convection independent 
of frequency exists, then the phase angle of the cross spectral density function is pro- 
portional to frequency. 

In general, if two stations are fairly close to one another and the disturbances 
at the two locations are caused by the same source, it is expected that their auto- 
correlation functions and auto power spectra are of the same general shape at one 
station as the other. A plausible model of the correlation function for any two Stations 
1, 2 was proposed in Reference 30 (Equation 5-8, p. 43) as: 


*12 ^ r ’ T ^ C 1 °2 ® 


-<J| 


V- ) 


(4-27) 


where a > a denote the rms values (with zero mean), r is the distance between the 
1 2 

two points, 6 is the spatial decay coefficient, t is time, and f 1 relates to the time of 
convection of the pressure from one point to the other. p q denotes the normalized 
auto-correlation function. Based on the above assumption of identical or near-identi- 
cal auto-correlation functions (different only in amplitude) , the phase angle 6 of the 
cross spectral density $ 12 and the cross-correlation coefficient as defined by 
Equation (4-14) may be expressed in the following: 


9 = -2nf t' 


(4-28) 



(r,r')=e 



y 12 ^ r ’ * ) 


(4-29) 


In Equation (4-29), y ^ is the square root of the coherence function defined by 
Equation (4-13). Apparently, is assumed to be independent of frequency f. Based 
on (4-29), the spatial decay constant 6 may be determined. 

Equation (4-28) shows that for a constant r 1 corresponding to a specific convec- 
tion speed of the disturbance, the phase angle is a linear function of frequency passing 
through the origin of the 9 -f plot. For this case, the convection speed = r/r 1 . In 
the following, correlation data based on pressure transducer pairs (1, 2), (2, 3), (2, 4), 
and (1, 6) as presented in Figure 63 to 66 are discussed. In these four figures, each 
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has four plots. The two plots on the left side give the modulus and phase angle of cross 

PSD 0 (f) as shown in Equation (4-10). The two plots on the right side give the co- 

xy 2 

herence function y (f) and the correlation function R (t) respectively. The cross 
PSD modulus has been normalized and the normalization factor is given on top of the 
plot. 


Referring to Figure 63, the distance between transducer pair (1, 2) is approxi- 
mately 22.86 cm (9 in.) along the chordwise direction. The shape of the cross-correla- 
tion function is similar to the shape of the auto-correlation function of either transducer. 
The phase angle 6 of cross spectrum density function is approximately proportional to 
the frequency with t 1 = 0.002 second. The t 1 data is consistent with the time delay based 
on the major peak in the cross-correlation plot. The above observations indicate that 
for transducer pair (1, 2), a convection velocity exists and the V value is approxi- 
mately 114.3 m/sec or 375 fps, with p = 0. 814. A similar conclusion may be drawn 

max 

for transducer pairs (2, 3) and (2, 4) located at the same semi-span station. Figure 


max 


= 0.710. The corre- 


64 shows t' = 0. 004 second for transducer pair (2, 3) with p 

sponding convection velocity is 57.2 m/sec (188 fps). Figure 65 shows t' = 0.01 

second, with r = 57.15 cm (22.5 in.) for transducer pair (2, 4). The corresponding 

convection velocity is 57.2 m/sec (188 fps), which is identical to that of transducer 

pair (2, 3). The p value is 0.567 for transducer pair (2, 4). 
max 

Between two pressure locations, the mean of the square root of the coherence 
function y x y(r) may be determined by an averaging process covering a certain fre- 
quency range. In case there is a dominating convection effect, then according to 

Equation (4-29), the mean of y (r, f) may be determined by the normalized cross- 
ly 

covariance function as defined in Equation (4-14). Based on this method, and using 
pressure data obtained in Run 2 of Flight 871, the contours of equal y are plotted in 
Figure 67 for Station Number 2. The contours are plotted only in the area where the 
convection effect is observed so that the approximate formulation of Equation (4-29) 
is applicable. The curves are left broken at locations where the convection effect 
is not observed. Figure 67 indicates that a distinct convection effect existed in the 
area near the wingtip. This effect is most prominent along the chordwise direction. 
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Figure 63. Correlation Between the Pressure Transducer Pair (1, 2), 
M = 0.925, h = 10,668 m, = 4°, 6 . = 12° 





smnaow 


S33y93U - 3SVH d 


00 

I 

<N 


aS 


<N 


2 00 

£ s 

c - 
0) o 

0) r-l 

I « 

0 J3 

pq ■= 


a> 

u 

u 

o 

U 


to 

« 

3 

bX3 


112 



Figure 65. Correlation Between the Pressure Transducer Pair (2-4), 
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Figure 66. Correlation Between the Pressure Transducer Pair (1-6) 
M = 0.925, h = 10,668 m, 6 = 4°, <5. - 12° 



Figure 67. Contour of Mean Square Root Values of Coherence Function y for 
Buffet Pressures Obtained in Run 2, Flight 871, xy 

M q = 0. 925, h = 10, 668 m, 6 = 4°, <5 f = 12° 
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4.6 Statistical Distribution 


The definitions of the probability density function p(x) and cumulative probability 
distribution function P(x) are given by Equations (4-16) and (4-17) respectively. An 
estimate of p(x) is obtained in digital form by dividing the range for x into an appropri- 
ate number of class intervals and by tabulating the number of observations and the per- 
centage of data in each interval. An estimate P(x) is obtained as the sum of the number 
of observations and the percentage of data from the class intervals. 

Typical statistical nondimensionalized probability density and distribution func- 
tions of pressure transducers 1, 6, and 23 are presented in Figures 68 to 70. The 
ordinate of a nondimensionalized probability density function plot represents the ratio 
of the event counts falling in a pressure level internal to the total number of event 
counts. The data are based on Flight 825, Run 7 low-frequency digital tape. The time 
span covers the period from 364. 5 to 366. 0 second when separation was fully developed 
over the upper wing surface. The time increment is 0.001 second. In general, the 
pressure data do not follow a normal distribution as would be the case for a purely 
random natural phenomenon. The data do not have a zero mean within the chosen time 
interval even though the complete batches of data representing the maneuver have zero 
mean. 
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SECTION 5 


ANALYTICAL PROGRAM 

The dynamic buffet pressure applied on the aircraft is quite complicated. Its 
origin may be traced to the turbulence in the separated and, possibly, unseparated 
flow regions, as well as to the fluctuating pressure characteristic of the unstable 
shock and the movement of the separation boundaries. In case the buffet pressure 
may be assumed to be stationary or near-stationary within a certain time period, and 
the auto- and cross-power spectra of the pressure exist,* then an analytical procedure 
may be applied to compute the aircraft response due to buffet. The analytical response 
data may then be compared with the flight test response data to examine the validity of 
the analytical procedure. 

The analytical procedure is applied in two different approaches. In the first 
approach, the data covering a substantial time span of the transonic maneuver are 
used which are assumed to be stationary. In order to correlate the average response 
behavior in the late phase of the maneuver including wing rock, the data covering the 
initial phase, i.e. , prior to buffet onset are excluded. This averaging approach is 
considered to be more suitable to low Mach number case (M Q = 0.75) where the large 
amplitude local shock oscillations do not exist. In the second approach, the time span 
covering the transonic maneuver is divided into a number of consecutive time seg- 
ments. Within each time segment, the fluctuating pressure data are assumed to be 
stationary. The aircraft response is then computed based on the same linear response 
transfer function and the segmentwise stationary buffet loads. 

The second approach mentioned above takes into consideration the ever changing 
buffet loads at various phases of the maneuver. It has the disadvantage that the time 


*An alternate condition is that the Fourier transforms of the pressure data at all 
collocation stations exist and the spatial correlations between any two stations are 
specified. Two extreme cases corresponding to the latter condition are described in 
this section and in Appendix I of the Report. 
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span of each segment becomes so small that the normalized standard error e of the 
power spectra is large based on a given equivalent resolution bandwidth B e (see 
Equation 4-22). For the aircraft under consideration, the more significant natural 
modes are located in a narrow frequency range. In order to resolve the closely packed 
modes, a small frequency increment 4f is necessary. Since B e = 2df in the applied 
processing technique (see p. 290, Reference 47), some compromise has to be reached in 
order to obtain a reasonable resolution and an acceptable level of spectral data error. 

In the following subsection, alternate methods of organizing the buffet pressure data 
and the results of varying bandwidth are described. The subsection is followed by 
a general description of the aircraft response transfer function, subsequent to which 
the response data obtained by the two analytical approaches are presented. 

5. 1 Force Input Based on Measured Buffet Pressure 

In order to compute the aircraft responses, force input based on the measured 
buffet pressure data is organized in matrix form using either the Fourier transform 
functions or the auto- and cross-power spectra. To assemble the force input data, 
the wing planform is divided into a number of subareas, with the pressure transducers 
located inside, as shown in Figure 71. The power spectral densities of pressure are 
assumed to be uniform over each subarea. Perfect spatial correlation is assumed 
between any two points lying within a subarea. In the following, the term "pressure" 
refers to the pressure over a subarea. 

Two methods, the Fourier Transform Method and the Correlation Method, are 
available to generate the pressure input data for the frequency response calculation 
program. Theoretically, the Correlation Method will yield the exact response solu- 
tion to the problem since it accounts for the actual spatial correlation of pressures. 

But the pressure input involves a 20 x 20 (total number of transducers on the top sur- 
face) Hermitian matrix whose elements are the auto- or cross-power spectral densi- 
ties of the corresponding pressure transducer pairs, for each given forcing 
frequency. As a result, tremendous amounts of data have to be generated and 
sorted. The Fourier Transform Method generates only 20 Fourier transforms of 
pressures for each forcing frequency. The auto power spectrum is formed by the 
complex conjugate product of the Fourier transform function. Based on the Fourier 
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transform functions of the subarea pressures, the corresponding power spectral 
matrix for the pressure input may be generated for two extreme cases. These are: 

1. Two subareas are perfectly correlated. For this case, the elements of 
the spectral matrix representing the pressure input are formed by the 
complex conjugate product of the Fourier transform functions of the two 
sets of the pressure data. 

2. Two subareas are uncorrelated spatially. For this case, the spectral 
matrix is a diagonal matrix whose diagonal elements are the auto power 
spectra of the subarea pressures. 

In reality, the spectral matrix of the buffet pressures should be somewhere 
between the two extreme cases. Because of the large amount of computation needed 
to deal with the spatial correlation of any two subareas individually (see Section 4), 
the analytical computation is limited to the extreme cases with which the experimental 
results are compared. 

In order to determine the normalized standard error and its effect on spectral 
data, PSD's are generated in two ways for Pressure Transducers 1, 4, 6, and 18 
data corresponding to Run 7, Flight 825 (M Q = 0.75, <5 n /<5f = 0°/0°). 

Firstly the PSD data are generated using a low frequency digital tape covering 
a time span T of 4. 096 seconds (362. 00 - 366. 096). The low frequency tape has a 
frame rate of 1000 per second, i. e. , the time between samples is . 001 second. A 
new set of data is then generated from the data samples using a 0. 002 time incre- 
ment. (The new time increment is an integer times the time between samples in 
order to obtain a desired frequency resolution.) The effective resolution frequency 
B e used in the process is 0. 9766 Hz. The normalized standard error is then 


f =. (B e T) 1/2 = 0. 500 

which is very high. The PSD's for the same data recorded on a higher frequency 
digital tape are also generated. The high frequency tape has a frame rate of 5000 
per second. The cut-off frequency of this tape is 2500 Hz. With T = 4. 915 seconds, 
B e = 9. 766 Hz, the normalized standard error is: 
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The PSD data for pressure stations 1, 4, 6, 18 obtained from the high and low 
frequency tapes are plotted one over the other in Figures 72-75. (Only the low 
frequency portions of the high frequency PSD are plotted.) Referring to the figures, 
except for lack of certain details in the high frequency PSD in the low'-frequency 
region, the high-frequency and low-frequency PSD's for any one set of pressure data 
are considered consistent. This observation gives credence to the use of the low 
frequency (high resolution) PSD's in the response study in case a certain amount of 
the normalized standard error can be tolerated. 
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Figure 71. Locations of Pressure Transducers, Accelerometers, and Control Points 
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Figure 74. Low Frequency PSD Plot of Pressure No. 6, 
M = 0.75, h = 7, 772m, 6 =0°,6 C =0° 
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5.2 Aircraft Response Function 


The aircraft is assumed in a basic configuration, i.e. , full fuselage tank, gear 
up, bare wing except that AEM-9B missiles are attached to both wingtips. In the 
symmetrical frequency response computation, a total of 40 control points are used 
which represent the normal deflections of the aircraft. The detailed control point 
geometry is given in Figure 76. In addition to those listed in this report, the weight, 
inertia, and stiffness data of the aircraft may be found in Reference 32. 

A total of 12 free-free flexible modes are utilized in the analysis. The highest 
eigen frequency is 54.19 Hz. Two rigid body modes are considered, i.e. , vertical 
translation and pitch rotation. The structural modal damping coefficient is assumed 
to be 0.05 for all the modes. The above modal data are used for response computation 
dealing with Run 7, Flight 825, M q =0.75 case. Simplified modal data are used for a 
nonstationary response computation dealing with Run 2, Flight 871, M q = 0.925 case, 
The latter case is described in subsection 5.4 of this report. 

For a number of frequencies within the band of interest the aircraft response 
functions are computed using a modal approach. For this purpose, the basic elastic 
modes of the aircraft are generated, and are used as the normal modes in the 
response formulation. 

The equations of motion for the determination of the frequency response functions 
are programmed using both the cantilever and the free-free modal solutions which are 
synthesized (see Appendix I). The frequency response functions are determined for 
deflection, acceleration, control point load, shear, bending moment, torque, and 
stress at specified locations. In this subsection, the basic building blocks used in the 
aircraft response analysis are described. 

Mass Stiffness Model 

The aircraft is represented by a system of control points with concentrated 
masses. Additional intermediate. masses are suspended on massless bars connecting 
the control points (Figure 76). The control point masses and intermediate masses are 
calculated such that the mass and inertia properties of the system are matched with 
those of the aircraft. The control point and intermediate masses are used to derive 
the coupled mass matrix for a dynamically equivalent system. Additional control 
points with no associated masses are added at locations which are required to correlate 
the response data with measured flight test data. The elastic axis of the wing is the 
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Figure 76. F-5A Wing, Fuselage, and Horizontal Tail Control Point 
Geometry for Symmetrical Analysis 



35% local chord near the root and moves toward the 40% local chord near the tip. The 
control points are located at 25% and 75% local chord locations as shown in Figures 
71 and 72. The flexibility matrix is calculated and merged for each of the three main 
sections of the model: right wing panel plus tip-store combination, forward fuselage 
and aft fuselage plus horizontal-tail combination. 

Modal Data for the F-5A Fighter 

The natural vibration mode data based on the mass-stiffness data described 
above are computed. (See Appendix n for computer program organization.) The fre- 
quency and mode description are tabulated in Table 5. Also listed in the table are 
typical ground vibration natural frequency data for the same aircraft configuration. 
The detailed flexible modal data are plotted in Figures 77 and 78. Only symmetrical 
modes are presented. 


TABLE 5 

THE F-5A FIGHTER NATURAL FREQUENCY MODE DATA 


MODES 

F(Hz) 

GROUND 

VIBRATION 

TEST 

FREQUENCY (Hz) 

DESCRIPTION OF MOTION 

1 

4.049 

4.45 

1st Wing Bending, Fuselage Bending (wingtip and 
fuselage nose are out of phase) 

2 

6.522 

6.60 

1st Wing Torsion 

3 

8.378 

10.2 

1st Fuselage Bending, Wing Bending (wingtip and 
fuselage nose are in phase) 

4 

16.983 

17.20 

1st Horizontal Stabilizer Bending 

5 

18.850 

18.40 

2nd Wing Bending, Fuselage Bending (wingtip 
and fuselage nose are out of phase) 

6 

21.989 


2nd Fuselage Bending, Wing Bending (wingtip 
and fuselage nose are in phase) , 

7 

31.037 


3rd Fuselage (Forward Fuselage) Bending 

8 

36.620 


2nd Wing Torsion, Forward Fuselage Bending 

9 

40.048 



10 

45.880 



11 

52.801 



12 

54.192 


. 




131 

















Figure 78. Higher Symmetrical Modes of F-5A 
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As described previously, the buffet pressure data were acquired on the right- 
hand wing. No comparable data were available on the left-hand wing. Furthermore, 
the wind-up turn did induce certain lateral motion of the aircraft. The extent of par- 
ticipation of the antisymmetrical structural modes in aircraft buffet response was not 
clear. No conclusive information was obtained from the left-hand wingtip accelerom- 
eter because of the variation in types of accelerometers used and the excessive noise 
picked by the left side accelerometer. It is due to lack of information, rather than 
theoretical considerations, that only the symmetrical modes are used for the aircraft 
response computations. 

Aerodynamics 

The oscillatory aerodynamic influence coefficients (AlC's) are computed based 
on unsteady aerodynamic theory assuming thin airfoil lifting surfaces and slender 
bodies. The AIC's relate the control point deflections whose time derivatives are the 
downwash velocities. Aerodynamic forces expressed in terms of AIC's for motion are 
given in Equation (1-8) of Appendix I. The resulting AIC's are functions of planform, 
Mach number, and reduced frequency. 

The AIC's for wing and tail motion are based on lifting surface theory. The 
basis of the lifting surface theory is found in the aerodynamic integral equation that 
relates the pressures on the surface to the downwash at the surface. A number of 
solutions to these integral equations by collocation methods have been proposed. The 
methods are based on certain approximate treatments of the kernel of the integral 
equation. In the subsonic range, the kernel of Watkins, Runyan, and Woolston is used 
(Reference 41). The method of solution is based on the expansion of the kernel function 
into a fifth-order power series in reduced frequency and was originated by Runyan and 
Woolston (References 42, 43, and 44). It consists of adjusting the'amplitudes of an 
assumed series of modes for the pressure corresponding to the downwash produced at 
the control points by a system of chordwise replacement loads. 

Tail downwash due to wing lift is accounted for (see Reference 45). The wing 
and tail AIC matrices are assumed to be decoupled. 

The AIC's for the fuselage due to motion are derived from the slender-body 
theory (References 43 and 46). The method employs the momentum theory of Munk 
and Jones (Reference 3), The disturbance flow component considered in the momentum 
computation is two-dimensional and is located in a plane normal to the direction of 
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flight. These fuselage AIC's are independent of Mach number. In the following para- 
graphs, the buffet pressures based on measured fluctuating pressure data are defined. 
The natural modal data and the unsteady aerodynamic coefficients are combined to 
form the aircraft transfer functions subject to the buffet pressures. 

In subsection 5.1, the spatial correlation and its effect to the pressure PSD 
matrix have been discussed. Another factor to be considered involves the separation 
of the measured pressure into proper components for analytical correlation purpose. 
Specifically, the measured dynamic pressure distribution may be separated into two 
components: the pressure induced by wing motion, and a dynamic pressure due to 
flow separation, etc., which can be considered as the genuine buffet pressure. The 
modal force components corresponding to the pressure components are indicated by 
T , f respectively. The measured force f m due to these pressure components can be 
written as : 




where 0 is a fraction factor defining the portion of induced pressure to be assigned to 
the subarea under consideration. For instance, if the induced pressure is assumed 
to be distributed evenly on the top and bottom wing surfaces, and if f m is composed 
only of upper surface pressures where the significant buffet pressure appears, then 
0 = 1 / 2 . 

The upper-surface force induced by wing motion is: 


d <f. / = 0w 




where j a | represents the modal amplitudes and |a 
T he matrix equation of motion is then: 

[z-(l- S )„ 2 A]{*}={fJ 


(5-2) 

is the generalized AIC matrix. 

(5-3) 


where Z is the modal impedance matrix and u> is the circular frequency. The detail 
mechanization of Equation (5-3) is described in the next subsection and in Appendix I. 
Substituting Equations (5-2), (5-3) into Equation (5-1), the genuine buffet force due to 
flow separation may be obtained as follows: 



5.3 Aircraft Response Analysis Considering a Major Portion 
of the Maneuver as Stationary 

The aircraft response functions are computed for a number of frequencies within 
the band of interest using the modal approach. The dynamic equations used in the 
response computation are matrix equations for a damped elastic system considering 
the buffet pressure input and the induced aerodynamic forces. Approximate solutions 
minimizing the errors for all degrees of freedom are obtained applying the Galerkin 
method. 

The computer program mechanization of the spectral response calculation is 
given in Appendix n. The computation starts with generating the symmetrical free- 
free modal data of the aircraft upon receiving proper idealized structural data as input. 
For a given Mach number, the aerodynamic influence coefficient matrices for wing, 
tail, and fuselage are generated separately corresponding to a sequence of reduced 
frequencies. The measured real-time pressure data are converted into Fourier trans- 
forms in the frequency domain. The frequency response program generates the 
Fourier components as well as the power spectral data of the response quantities based 
on the available modal, AIC, and pressure data. Assume that the measured pressures 
are random and ergodic. Two extreme cases of the spatial correlation between pres- 
sures are considered in the response computation as described below. 

1. Perfect spatial correlation condition: For this case the coherence functions 
for pressures always equal unity. The frequency response computer pro- 
gram computes the Fourier Components of the response quantities. The 
corresponding PSD values may be obtained as follows: 

I 2 

F r W | » (5-5) 

where T is the time span and F r is the Fourier component of the response 

quantity. Equation (1-27) in Appendix I is a special case of Equation (5-5) 
where F r (w) = (-u /G) h^(w) is the Fourier transform of acceleration at 
control point i. 

The disturbance pressures due to flow separation may be computed based 
on Equation (5-4). 
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2. 


Spatially uncorrelated condition. For this case the modal force PSD 
matrix |"^ f l is formed by a diagonal pressure PSD matrix. The coherence 


functions of the pressures in two subareas are always equal to zero. 
PSD's of the response functions are generated directly. The matrix 
equation equivalent to Equation (5-3) is of the following form: 



z-(i 


-/3)cu 2 A 


rl 


<t> 


:] [jz*-(l-/3)cu 2 a] ] 


The 


(5-6) 


H(V) 




In Equation (5-6), <f> (u) is the matrix of the modal response PSD. Equation 
(1-28) in Appendix I is identical to Equation (5-6), except that a factor (u /G) 
is added to convert the displacement PSD matrix into an acceleration PSD 
matrix. 


In case (2), the measured force cannot be physically separated into two compon- 
ents corresponding to Equation (5-1), as for case (1) where the pressure data were 
perfectly correlated in a spatial sense. For the actual case where spatial correlation 
between pressure subareas exists to a certain degree depending on the flight condition 
and the relative locations of the pressure subareas, Equation (5-6) still applies. The 
exception is that the matrix <^J for the latter case is formed by a full-pressure PSD 
matrix whose elements are the Fourier transforms of the auto- or cross-correlation 
functions. 


Using the above formulation, the analytical computation of aircraft response 
functions was performed for the following configuration: 


Flight number 
Mach number 
Altitude 
Flap settings 
Total time span 

Low frequency digital tape frame rate 
Time increment 
Frequency increment 
Spectral frequency range 


825, Run 7 
0.75 

7,772 m (25,500 ft) 
( 0 °/ 0 °) 

362-367 sec 
0. 001 sec 
0. 002 sec 
0.488 Hz 
3.5 - 46.5 Hz 
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The comparison between the flight test response data and the computed response 
data of the aircraft are based on the PSD of acceleration at key locations including the 
forward wingtip, the aft wingtip, and the C,G ti Figures 79-81. In these figures the test 
data are represented by heavy solid lines. The analytical data based on the perfectly 
correlated pressure are shown as dotted lines, while those based on spatially uncor- 
related pressures are represented by light solid lines. The differences in PSD's 
between the two sets of analytical data are indicated by the shaded area. 

Referring to Figures 79, 80, the differences of computed PSD's based on zero 
correlation and perfect correlation are in most cases limited to about a factor of 10 db. 
In general, the computed data based on zero correlation are somewhat closer to the 
PSD based on flight test response data. A major deviation between test and computed 
response PSD's appears in the frequency range above 25 Hz. In this range, the test 
data indicate substantially higher bending type vibrations, while the computed data in- 
dicate high amplitude torsional vibrations. In Figure 81, the substantial C.G. responses 
to the primary wing bending and torsion modes (at 4. 049 and 6. 522 Hz) are not properly 
reflected in the analytical data. 

The large deviations between the analytical and test response PSD's are attributed 
to a number of factors. These are listed below: 

1. The transonic maneuver is a nonstationary process. The 5-second interval 
used in the present computation does not take into consideration the non- 
stationary dynamic effect apparent during the interval. This factor is reme- 
died by the segmentwise stationary computation carried out in a subsequent 
subsection. 

2. The analytical data are computed based on the assumption that the buffet 
pressures on both wings are symmetric about the center section of the 
fuselage. The nonsymmetrical components of the actual buffet pressure 
distribution, if any, are ignored. The peaks corresponding to antisymmet- 
rical modes are thus missing from the analytical plots. Furthermore, the 
dynamic contribution by the tail surfaces can be significant and is ignored. 

3. The spatial correlations of pressures are neither perfectly correlated nor 
noncorrelated. The degree of spatial correlation is dependent on the trans- 
ducer locations. For example, the measured pressures near tip locations 
are highly correlated, while the pressures in the other region, especially 
those near the wing root, seem to be lightly correlated or uncorrelated. 



WINGTIP (AFT) ACCELERATION PSD - G'VHz 


TEST DATA 


COMPUTED DATA 

ZERO CORRELATION 

COMPUTED DATA 

PERFECT CORRELATION 
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I (VI 


8 9 10 


FREQUENCY - Hz 

Figure 79. Experimental and Analytical Data of Right-Hand Wingtip (Aft) Acceleration 
Power Spectral Density, M q = 0.75, h = 7,772 m, $ n = 0°, <5 f = 0°, 

Real Time Data Based on Time Span (362-367 sec) 




WINGTIP (FORWARD) ACCELERATION PSD -G 2 /Hz 



FREQUENCY -Hz 

Figure 80. Experimental and Analytical Data of Right-Hand Wingtip (Forward) 
Acceleration Power Spectral Density, M q = 0.75, h = 7,772 m, 4^ = 0°, 

4^ = 0°, Real Time Data Based on Time Span (362-367 sec) 
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Figure 81. Experimental and Analytical Data of eg Acceleration Power Spectral 
Density, M q = 0. 75, h = 7, 772 m, <5^ - 0°, = 0°, Real Time Data Based 

on Time Span (362-367 sec) 
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4. The computed response data, especially the frequency shift at the resonance . 
frequencies, depend on the induced pressure contribution parameter of /? 
which is arbitrarily assumed. The frequency shift is sensitive to the fi 
changes because the imaginary part of the unsteady aerodynamic force is of 
the same order as the structural damping force which determines the response 
peak frequency shift. The real part of the unsteady aerodynamic force is 
usually substantially smaller than the spring force of the structure. Also, 
improvement can be made on the analytical AIC matrix. 

5. In the neighborhood of a natural frequency, the computed response PSD level 
depends strongly on the modal damping coefficient input gp, which is assumed 
to be 0.05 in the present investigation. Furthermore, in order to streamline 
the computation, simplified structural models are used. 

6. Normalized standard errors in PSD computation are introduced in both test 
and analytical data. 


Induced Aerodynamic Pressure in Buffet 

A subject of interest in a buffet test program is the magnitude of the induced 
aerodynamic forces (as a function of frequency) relative to the overall dynamic 
pressure force during buffet. In these subsections, the computed induced aero- 
dynamic pressure data based on the measured buffet pressure distribution and the 
aircraft responses are presented. 


For the case of pressure distribution with perfect spatial correlation, it is pos- 
sible to separate the genuine pressure p^ due to flow separation from the measured 
pressure p . Since the Fourier transforms of the total pressure p. induced by wing 
motion in a uniform flow may be computed, the relative magnitude of the induced 
pressure at the top surface of the wing to the measured oscillatory pressure p m may 
be determined by: 


= 0 

P, 

/ 

P™ 


1 


m 


1 - V p 


m 


(5-7) 


at different transducer locations as functions of frequency. Typical pressure ratio 
plots at Transducer Locations 1, 4, 18, 23 are shown in Figure 82. The peak loca- 
tions in the plots correspond either to a resonance frequency which induces a large 
deflection at that transducer location, or a very low value of the measured pressure. 
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In general, the ratio f is less than 0. 3. In other words, the genuine buffet pressure 
P b dominates in most frequency ranges, while the induced pressure is less significant 
in the response computation. The top plot of Figure 82 shows the Z values for 
Transducer No. 1 where the relative amplitude of the induced aerodynamic pressure 
yields a very high peak: 


f = 7. 324 Hz 
f = 1.086 

p m = (-7.31 + i 5. 81)N/M 2 , 
/?p. = (-2.12 - i 9. 92)N/M 2 , 
p b = (-5.18 + i 15. 73)N/M 2 , 



= 9. 31 N/M 2 
= 10. 14 N/M 2 



16. 56 N/M 2 


The highest peak occurs near the second eigen frequency where Transducer Location 1 
has a large deflection. For this case, even though £ is greater than unity, | p b | 
is still 1. 6 times the value of /3 1 j while the two complex pressure vectors /3p^ 
and p b are out of phase. 


5. 4 Aircraft Response Analysis Considering the Maneuver as 
Segmentwise Stationary 

In order to improve the analytical correlation of aircraft response results, the 
major portion of a transonic maneuver is divided into a number of segments. In each 
segment, the buffet pressure data are assumed to be stationary. The aircraft is then 
subjected to the consecutive application of the buffet loads and the cumulative dynamic 
effects are reflected in the time varying PSD of response. The detail formulation of 
the segmentwise stationary approach and the computation flow diagram are given in 
Appendix III of the report. 
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The approach is applied to the following set of flight test data: 

Flight number 871, Run 2 

Mach no. 0.925 

Altitude 10, 668 m (35,000 ft) 

Flap settings (4°/12°) 

Low frequency digital tape 1000 per sec 

frame rate 



Time increment 


0. 002 sec 


Frequency increment 0. 488 Hz 

Spectral frequency range 1. 4 — 20. 0 Hz 

The real-time data of this run are given in Figures 19-28. Altogether, 
dynamic pressure data covering the time span 073. 00-082. 10 second are used. 

The data, assumed to be perfectly correlated, are divided into four (4) equal time seg- 
ments. The computation is carried out using the rigid body plunging mode and the first 
three symmetrical flexible modes (see Table 5 in subsection 5. 2). The pitch mode is 
not included because of lack of tail surface dynamic loads data. Corresponding to the 
above flight condition, and considering plunge-pitch coupling, the natural frequency 
of the short period mode of the aircraft is 0.3947 Hz. The damping ratio of the 
short period mode is 0. 245. 

Figure 83 gives the time varying influence function of the aircraft corresponding 
to the first flexible mode (first wing bending) and during the first time segment. The 
absolute values of the influence function | 1^ ^(tjCj) j are plotted in G-level per unit 
modal force. (See Appendix III, Equations III— 12 and III-23, for definition.) 
Corresponding to a rectangular shaped time function e^t), the influence function is ' 
independent of frequency at the initial time. It peaks at the natural mode frequency 
as time passes. The function value drops rapidly in frequency ranges other than the 
natural frequency at the end of the time segment. This behavior illustrates the 
dynamic carry-over effect of the modal response in the subsequent time segment. 

Figure 84 gives the mean square values of the rigid body (plunging) mode and 
three flexible mode forces corresponding to the four time segments of the transonic 
maneuver. Only the mean square contributions within the frequency range 1.4-20.0 
Hz are taken into account. The relatively high modal forces in the second time seg- 
ment reflect the shock oscillations at buffet onset (see, for instance, Figures 19, 22). 
The modal forces diminish in values in the third time segment and reach their peaks in 
the last and fourth time segment when the angle of attack reaches a maximum. 

Based on the preliminary data such as shown in Figure 83, 84, the nonstationary 
response PSD's for the C.G. and two stations of the right wing tip are computed at 
specific times, one each within the four time segments. These results are plotted in 
Figures 85-87. Also plotted are the corresponding segmentwise stationary PSD's 
based on the flight test response data. For the two sets of response PSD's at the right- 
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hand wing tip, the computed data are too high in the first segment. This is believed 
due to the higher damping (of the Coulomb type) at the initial phase of the maneuver 

not accounted for in the computation. The correlation is more satisfactory in the 
third and the fourth segments when wing rock takes place. For the C.G. acceleration 
PSD's, Figure 87, the computed first segment PSD values are low because more 
spectral energy has been contributed to wing vibrations at this time segment as 
explained above. For the later segments, with some exceptions, a somewhat better 
correlation between the computed and flight test PSD data is realized. 

In Figure 88, the time varying mean square acceleration for the C.G. and the 
two stations of the right-hand wing tip based on the segmentwise stationary formula- 
tion are plotted. Again, only spectral contributions in the frequency range 1.4-20.0 
Hz are accounted for. The figure reflects the very low C.G. response in the first 
time segment as described above. The substantial differences in mean square accel- 
erations between the two wingtip stations reveal substantial wing torsion mode partici- 
pation. In general, the correlation of analytical and test response PSD data is more 
satisfactory using the segmentwise stationary approach as compared to the approach 
where a major portion of the transonic maneuver is considered stationary. 
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Figure 88. Time Varying Accelerations of the CG and Right-Hand Wingtip During the Transonic 
Maneuver, Run 2, Flight 871, M q = 0.925, h = 10,668 m, <5^ = 4°, = 12° 

(Analytical Data) 


Section 6 


CONCLUSIONS AND RECOMMENDATIONS 

In the previous sections of this report, the findings of the buffet flight test 
program, the processed data, and the analytical correlation results have been pre- 
sented. Typical buffet onset data of F-5A which may be found in previous literature 
(Reference 26) are not included in this report. 

The emphasis of the present program is on the acquisition of the dynamic buffet 
pressure and response data. Substantial effort has been dedicated to the spectral and 
statistical processing. The results have been presented in the report. Certain con- 
clusions and observations relating to the flight test and processed results have been 
given, together with the presented data. In this section, major conclusions and 
recommendations are outlined. 


6.1 Conclusions 

1. Precise and detailed dynamic pressure data acquisition for aircraft during 
transonic buffet is feasible and can be productive. No significant adverse 
effect of the added instrumentation and change of airfoil geometry to the 
natural flow pattern was observed based on the dynamic flow development 

on the wing and the comparison with the tuft data. Fairly consistent mapping 
of the separation boundaries was achieved. 

2. Shock instability in a steady-state uniform flow has been observed in 
previous wind tunnel and theoretical treatments under specific conditions 
depending on Reynolds number and local Mach number. The dynamic 
pressure of the F-5A showed definite oscillations where the separation 
boundary was located . It is uncertain whether the cause of the shock 
oscillation was solely due to the inherent shock instability, the cause might 
be partially attributable to the interaction with aircraft motion and angle of 
attack changes. This point remains to be resolved. In general, the 
amplitude and pattern of the shock oscillations were such that they were 
expected to contribute significantly to the aircraft motion. 
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The pressure PSD distributions of a buffeting aircraft were such that the 
establishment of a mathematical model based on a number of flight con- 
dition parameters was feasible. Whether the mathematical model is general 
enough to cover various types or designs of aircraft is a matter subject to 
future investigation. 

4. The spatial correlation of the buffet pressures at various locations of the 
wing was either high or moderately high in the outboard stations near the 
wingtip. Substantially lower spatial correlation was observed for pressures 
between inboard stations or between one inboard and one outboard station. 

5. Given detail buffet pressure data, the aircraft responses were predicted 
in the low-frequency range using linear transfer function technique . The 
segmentwise stationary approach was the preferred approach. With addi- 
tional refinements in both buffet pressure model formulation and aircraft 
transfer function computation, a practical method can be developed for 
aircraft buffet response prediction. 



6. 2 Recommendations 


Based on the observations and conclusions presented above, a number of 
recommendations are listed below which may be used as guidelines in future research 
and development work on aircraft buffet. 

1. It is most desirable to acquire steady-state pressure data for the same 
aircraft under identical flight conditions in order to reach a better under- 
standing of the shock and separated flow interaction as well as other aero- 
dynamic effects of a buffeting aircraft. 

2. The flight test and processed data presented in the report may serve as an 
excellent set of baseline data to study the Reynolds number and scaling 
effects in transonic wind tunnel testing where a scale model is used. Items 

» such as the effect of angle of attack change during data acquisition, the 

adequacy of using transition strips in low Reynolds number tunnel testing, 
and the method of correlation of scale model and full-scale data are of utmost 
interest and should contribute to the state of art in transonic wind tunnel 
test techniques. 

3. The extent of participation of the rigid body modes (symmetrical and anti- 
symmetrical) as well as the participation of the anti-symmetrical struc- 
tural modes in aircraft buffet should be investigated. Dynamic loads on 
the tail surface which are subject to wake flow are to be included. Analyti- 
cal and experimental work should be carried out simultaneously on this 
subject. 

4. It is considered fruitful to further develop a mathematical model of the 
buffet pressures for a type or family of aircraft. The spatial correlation 
of the buffet pressures is a subject worthy of additional work. Definite 
progress may be achieved in establishing and improving the analytical 
prediction techniques for aircraft buffet responses. 

5. Design refinements and/or flight maneuver techniques should be looked 
into with a view to delaying the buffet onset or to reducing the amplitude of 
wing rock during a transonic maneuver. 
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APPENDIX I 


METHOD OF FREQUENCY RESPONSE CALCULATION 

In this appendix, a detailed derivation of the equations of motion and the response 
quantities (deflection, acceleration) are represented. The equations are developed by 
the collocation method extended to the approximate modal solution using Galerkin's 
method. Using the modal solution, the frequency response functions due to an external 
sinusoidal force or the Fourier transform of a stationary ergodic random pressure 
field are determined. In the analysis, the following sign conventions are used; The 
coordinates (x, y, z) and the forces and deflections are positive forward, starboard, 
and downward. # 


The analytical computation of the response data starts with the determination 
of the natural frequency modes of the aircraft as a free-free body. The reduced 
eigen equation for a free-free system may be expressed in the following form , as 
shown in Reference 50: 

[k u + k u T M ;l T T M u ] J*,j= 4 K] K 1 4* 0 <M> 

where is the eigen vector matrix of the control points, is the stiffness 
matrix for the system with the reference point fixed (i.e. , no rigid body motion), 
is the mass matrix of the control points, T is the rigid body transformation matrix, 
and is the associated circular eigen frequency. The corresponding reference point 
mode shape x~ may be expressed as a function of X^ 


= ' M 22 TTM 11_ 


1*1 


\ 


( 1 - 2 ) 


where is the reference point mass matrix and is assumed to be a null matrix. 

Let and be the modal matrices whose columns are composed of the eigen- 
vectors x ^ and x q respectively. The modal matrices satisfy the following orthogonal 
conditions: 


k 1 


t m u Xi + xJm 22 x 2 


= fM F J 


a-3) 
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where is the generalized mass matrix in diagonal form . In the frequency domain 
the equation of motion of a damped elastic system free in space may be written as: 


-" 2 [ m u] l h l I + (‘“[c]+ [*„]) K -Th 2 j = jq <«>} <I-5> 

-J [t T M u ] I h l| -" 2 [ M 2 2 ] Kj - [ T ] T | f 1 <“>! + ! f 2<"») ff- 6 * 


where h and f represent the Fourier transforms of the deflections and the external 
forces respectively. The subscript 1 in the above equations refers to the control 
points while the subscript 2 refers to the reference point. C is the damping matrix 
for the system. Equation (1-5) is the equation of motion of the control points if the 
reference point is held stationary. Equation (1-6) is the equation of motion of the 
reference point when all inertia forces and external forces are accounted for through 
the use of the transformation matrix T. The external forces are assumed arising 
from the motion and from some external sources such as the buffet pressure, etc. , 
‘thus: 


f (<j) 


= i*i<«> 


+ f 




(1-7) 


In case the pressure distribution is measured through pressure transducers at 
various locations, the relationship between the measured pressures and the buffet 
pressures is defined in Equation (5—1) . Then Equation (1-7) may be expressed in an 
alternate form as follows: 

jf( W )} = (1-/9) |f. (<4} + |f m Mf (I ~ 7a > 

The forces due to the motion may be expressed in terms of the oscillatory AIC's: 



2,2 
pw b r 




where 


P = air density, 


Kg/m 3 



= reference semichord m 
= semispan length, m 
= aerodynamic weighting matrix 
= oscillatory AIC matrix 


( 1 - 8 ) 
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Based on the modal data obtained previously, the structure deformation may be 
represented as: 


ki -Wki + HI 


a 


R 


*MK 


+ 


“r i 


(1-9) 

a-io) 


where (jp and q>^ are the free-free and the rigid body modal amplitude respectively. 

Premultiplying Equation (1-5) by T T and subtracting from Equation (1-6) yields 
the equation of motion of the rigid components: 

-“ 2 [ s 22 ] kl - [ t ] t (•“ [ c ] + [ k n] )! h i - Th 2 j ={ f 2m <">| p- 11 ' 

Substituting Equations (I-7a) ■— (1-10) into Equations (1-6), (1-5), and (1-11) yields: 

! 1 ) 

(tJ )j + ( f 2m (w) ( (1-12) 


: 4 Jkl k 

-“ 2 [* n x i] + ( ‘“[ c ] + [ k i i] )[ x i - tx 2] 


~ j f lm (cu) ( + l e F 


a F r w 


M 11 T 


r W 


-(4 


-i r iT 


M 22 X 2 


( 

(i Jc]+ IkJ) 


(1-13) 


X 1 ' TX 2 


)ki 


U) 


M 


22 


KI* { ' 2m<“>! + ( £ 


R 


where 


[Kf] - [j M 11 X 1 + M 22 X 2 _ 

[ s rr] - [ T 


IT T M U T * m 22 


M = M + (l-£) pb‘ S 


w 


c h 


(1-14) 

(1-15) 

fl-16) 

(1-17) 


e„ | = errors 


R 


Applying the Galerkin principle to the equations of motion, (1-13) and (1-14), 
the following equation is established: 


[ x i] T kf + h] T l 


v* 0 


(1-18) 
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where 


Equation (1-18) leads to the generalized equation of motion: 

T 

■S 

P*ff] = 


f lm^ + 


= I M n X i + X 2 M 22 X 2 


(1-19) 

( 1 - 20 ) 


In order to simplify the problem the generalized damping matrix is 
assumed to be a diagonal matrix: 

T 


X 1 " TX 2 


If. f 
W F k Fj g Fj 


X 1 " TX 2 

k j b j » ^ L j , v , 

where gj^ is the structural damping coefficient in the £ tn free-free mode and k^, is 
defined in Equation (1-4). 

The rigid modal amplitude is found by solving Equation (1-12): 


(1-21) 


“r ”w 2 


M 


RR 


'f 2 


MrF 


“F + 


T ' f lm! + l f 2m 




(1-22) 


Substituting Equation (1-22) into Equation (1-19) leads to the final response equation: 


“F = 


{f i m r 


where: 


1[ B | 

= (^FF "^RF ^RR M Rf) + lW [ C Fj + f k Fj 


B 



T _ 

1 



— T — -1 T 
M M T 
RF RR 


B 


2 


X 0 

T _ 

— T — 1 

2 


RF RR 


(1-23) 

(1-24) 

(1-25) 

(1-26) 


The Fourier transforms of the control point deflections are then obtained from 
Equation (1-9). 

In case the pressure transducer locations are different from the control points. 
Equation (1-25) may be expressed in an alternate form as follows: 


1 = To 1 T _ 

r— t — -i ti 

MMO 

J [^FJ 

"RF "RR ''R 


where: 


q t 


and Q 


R 


are the flexible and rigid modal matrices respectively, 


corresponding to the pressure transducer locations. 
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For an ergodic random pressure field , the PSD of the control point accelerations 
depends on the spatial distribution of the pressures. In the calculation, two 
extreme cases are considered: 


1. 


Perfect spatially correlated pressure field: 



(1-27) 


where 4>r (u) is the PSD of acceleration at a control point, subscript i iden- 
n i 

tifies the control point, T is the time span, G is the gravitational acceleration. 


2. Spatially uncorrelated pressure field and f„ = 0. The condition f_ =0 

zm Ztcl 

indicates that buffet pressures are applied only at the control points and not 
at the reference point. 


where 


and 


\ (to) J = G2 r 


H(w) 0 


Mi 

lm 


* 

H («) 


(«>f <">)« 

H(w) 


TT I'lm/^l 2 a « 


= X, - X., M, 


1 ^2 RR RF 


A 


-1 




0J‘ 


X 2 M RR X 2 


(1-28) 

(1-29) 

(1-30) 


Equation (1-28) is an expression for the acceleration PSD. On the right-hand side, by 

4 2 

dropping the factor oj /G , the remaining portion gives the response (deformation) 
PSD of the aircraft. This latter expression is identical to the right hand of Equation 
(5-6). 
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APPENDIX E 


COMPUTER PROGRAM MECHANIZATION OF 
FREQUENCY RESPONSE CALCULATION 


The flow diagram to generate the free -free frequency response data of an 
aircraft is shown in Figure 102 . There are nine functionally independent programs 
involved as described below: 

1. Program F019 

It generates the flexibility matrices for wing, fuselage, and tail with 
reference point(s) fixed. Only beam elements with torsion and one-way 
bending stiffness are considered. The flexibility matrices are punched 
out by the program. 

2. Program RF040 

It performs the free -free modal analysis for the complete aircraft upon 
receiving as input the flexibility, mass, and rigid body motion matrices 
of all components. The program generates in punched-out data form the 
free-free mode shapes and the natural frequencies in Hz . 

3. Program AA 

It computes the corresponding values of reduced velocity I k 

of the frequencies of interest which are consistent with those of the pressure 
input generated by program SPECTRAL. Additional input is the aircraft 
velocity, the reference semi-chord, and the prearranged indices for the 
frequency input. 

4. Program RHM15 

It formulates the AIC matrices for the fuselage based on the slender body 
theory. Each computed AIC matrix which corresponds to a given reduced 
velocity / is stored on disc file Number 16 identified by its pre- 
arranged index number. 
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5. Program RF022 

It formulates the AIC matrices for the wing and the tail based on the 
unsteady lifting surface theory. The computed AIC matrices which 
correspond to a given reduced velocity ^k^j and a specified Mach number 
are stored on Disc File Number 15 for the wing and on Disc File Number 17 
for the tail , identified by the prearranged index number for that frequency 
|or k^j respectively. 

6. Program M0FIT 

Since the pressure transducer locations are different from the control 
points on the wing, the computer program "M0FIT" is developed to compute 
the corresponding modal displacements at the pressure transducer loca- 
tions for each flexible and rigid-body mode. Corresponding to each mode, 
the program defines the deformation curve of the elastic axis of the wing 
based on the least square of a polynomial up to the third order: 

h = A 0 + A x y + A 2 y 2 + A 3 y 3 

The chordwise movements are assumed to be a rigid body rotation about 
the elastic axis. 

Program M0FIT punches out the output matrices [aJ[Q f J and [Aj [Q^J , 
where [Qp] and [Qp] are the flexible and rigid-body modal matrices 
corresponding to the transducer locations respectively, and TaJ is a diagonal 
matrix whose elements represent the wing subarea corresponding to the 
transducer locations. The products of matrices fAj [Q f ] , f AJ [Qp] and 
the Fourier transforms of the measured pressures yield the appropriate 
modal forces acting on the wing. 

7. Program R009H 

It calculates the frequency response quantities (deflections, stress, etc.) 
due to external sinusoidal forces (i.e,, the Fourier transforms of the external 
forces at that frequency) based on the Galerkin's method. The general input 
includes the free-free mode shapes, the corresponding natural frequencies, 
the rigid body modes, the modal damping coefficients, the matrices FaJ 
[Qp] and [ Aj [Qp] from program M0FIT, and some parameters such as 
p, V, b r> the reduced velocity, the corresponding S, etc. The frequency- 
dependent input includes Fourier transforms of the measured pressures at 
transducer locations and the AIC matrices which are read directly from 
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the disc files. Two kinds of pressure spatial correlation are considered: 
perfect correlation and zero correlation. 

8. Program SPECTRAL 

The spectral analysis program generates the Fourier transforms and the 
• power spectral data for all the recorded pressures and responses. The 
program receives the low-frequency digitized tape as input. It plots the 
output data listed below: 

a. Fourier transforms 

b. Auto- and cross-power spectral densities 

c. Auto- and cross-correlation functions 

d. Coherence functions 

The Fourier transforms of the pressures are punched out in a consecutive 
order with a specified frequency increment. 

9. Program FTPN 

The program reads in and sorts the punched output (Fourier transforms 
of pressures) from program SPECTRAL. Its output is another set of 
punched cards in a specified format consistent with the input requirements 
of the frequency response computer program R009H. 
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APPENDIX m 


MULTI-MODAL RESPONSE OF AIRCRAFT TO SEGMENTED 
STATIONARY RANDOM EXCITATION 

The formulation is an expansion of a nonstationary response analysis developed 
in References 51-53. The origin of the references may be traced to still earlier work 
on nonstationary analysis (for instance, Reference 54). The input function is assumed 
to be the sum of n functions defined in consecutive time segments. Each function is 
the product of e r (t) and f^ r (t). e r (t) is a deterministic function. In the present study, 
it is assumed to be unity within the time span (t , t^) and zero elsewhere. f, r (t) 
is a stationary random function. The development of the multi-modal response due 
to the segmentwise stationary input modal forces Q^(t) are described below. 

I. Modal Forces 
n 

Q/t) = e r (t) f fr (t) (m-1) 

r= 1 

i = mode index 

Cross correlation of modal forces Q, and Q 

t m 


%Q k (t l’ V 

^>k 

r,s 


n n 

S e r (t l ) e s ( v R f £r f ks (t l “ V 
r=l s- 1 ? i 

mode indices 

time segment indices 


(in-2) 
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Cross PSD of modal forces and Q^: 




(u v u 2 ) = -ig ff R Q^ Q k (t l’V exp [ ' u 2 t 2 ) ] dt l dt 2 

4 n JJ- oo 


n n 


El/' 


r= 


1 s=l *'-« 


S, „ (w) A (cj - <j ) A* ((j - u 0 ) dw (III -3) 

i . I, r j. s <s 

fr ks 


w l ,a> 2 = fre< l uenc y coor di na tes 


A r (« - w i) = ^ / e r (tj) exp £i (u - w x ) t ± 


dt-, 


exp [i (w - w 1 )t rb ] - exp [i(« - w 1 )t ra ] 
2-rri (w - w^) 


A (w - u,) is the Fourier transform of e (t..). 
r 1 r l 

n. Segmented Stationary Modal Forces 

Modal forces for r ^ segment 

{ f r (t)} = [xf Taj |p r (t)| 


where 


(ni-4) 


(m-5) 


[x] = [[Cypy], modal matrix corresponding to pressure transducer locations 

[Q f ] = free -free mode shapes 

[Q r ] - rigid body mode shapes 
|p(t) j = measured pressures 
l~Aj ~ subareas associated to transducers 
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PSD matrix of segmented stationary modal forces f r and f g 


s. . (u)l - [xfrAj[s (")1 T A J [x] 
L Vs -1 L J 


(m-6) 


where 



‘ 1T r) P P s ,<“- aT S >' 


n 



perfect spatial correlation 
zero spatial correlation 


i.j 

4T 


F p(^, 4T) 


indices of spatial locations 
time span for each segment 

Fourier transforms of measured pressures over a finite 
period 4T. 


(m-7) 


HI. Modal Responses 

PSD matrix of modal coordinates a 


[s„ (« 1 . U2 >] - [H( Ul >] [S Q ( Ul , <.,,)] [H«(- 2 )] T 

n n 

- ZZ [«<"!)] W T rA J [\p s <«>] TAj 

r= 1 s= 1 -oo 

[X] [H*(w) 2 ] T A* (w - w 2 ) dw (HI- 8) 

where 

[H(w)] = the modal transfer function matrix. Similar matrix is 
shown in Equation (5-6) in the stationary analysis. 

IV. Structural Response 


Structural deflection response at a specified location (i.e., accelerometer location) 

w = lyj { Q } . (ni-9) 

where 


Lyj = 


the modal shape matrix corresponding to the specified location 
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Power spectral density of w and mean square value w (t) . 


n n 

- » x ) |Y| [H(^j)] [X] T TAJ [s^ ( U )j 

r= 1 s= 1 - oo 


fAJ [X] [H*'(o> 2 )] T {y T }a*(o> - u 2 ) do> 


(in-io) 


w 2 (t) 



S (u 
w' 1 


w 2 ) exp 




dw^ dt «> 2 


"SZ /“L y j[‘r (t -“> 

r=l s=l -oo 



where 





] T M 


dtj 


J°° (H( Ul )j A r (w - w x ) exp (i^t) d Wl 


(m-ii) 


(in-12) 


an influence Unction corresponding to time segment r. 
elements of the modal matrix [H(u)] as shown in section V below 

i,j = row and column indices 
V. Aircraft Modal Transfer Functions 

* 

The dynamic equation of the aircraft in terms of the free-free and rigid body 
modes is: 



- 2 [ M FF] + r ^F^F + ig F W )J 

-“ 2 [ m rf 

] r 


;" 2 [“rp] 

1 

w 

l 1 


|{ a R|j 


|[Q f ] T f A jl F p[ 

[Qr] T fAjjFp} 


(m-i 3 ) 
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The modal transfer function matrix is: 




-" 2 [ m ff] + rw F + 


-“ 2 

-" 2 [“re] 


(m-i4) 


where 


I°r1 

I F P} 


free-free modal amplitudes 

rigid body modal amplitudes 

Fourier transforms of the measured pressures 
generalized masses 

non-zero natural circular frequencies 


\ 


g = modal damping coefficients 

r 

For [M 22 ] = 0, [M 22 J = 0 case, J^T T xj = 0, the modal mass terms are 
defined below: 


[ m ff] 


X I M U 



[ m rf] " 
[ M Rr] 


T T M n xj 

t t m 1x t] 



AIC included 


(ID- 15) 


The detail derivation and the definitions of the terminologies used in Section V 
are given in Appendix I of this report related to stationary response analysis. 
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VI. Modified Impedance functions 


Assuming that no aeroelastic coupling exists between the rigid body modes and 
the free-free modes, the impedance function matrix may be rewritten as follows: 




-1 


-J r M ff j + + ig pW )j 0 


Z r(-) 


(m-i6) 


In actual computation the following expressions are used: 
[M FF <*>] - “ <1 V, ( g A <"/»), 


(HI-1 7) 


1 +( t + ifi t 

It is noted that the modal aero-damping coefficient at w = approximately equals to 
(<^ + i/r )/(l - p). C £ and ^ will be used later in Equation (III-19). 

Z R (w) is a (2 x 2) matrix corresponding to the downward displacement and the pitch 
angle. The elements of Z R (w) are: 


! R (">| 

[ Z R.<“ ) | 




V"> 


Z rH 


12 


21 


22 


- u c / 4m _ r \ io> r 
2 y2 \pSc ” zb/ + V zb 

1 TT 2 0 / u 2 c „ iw _ \ 

P v Sc I o c • + TT C I 

2 \ 2 y2 ma V mo J 


! o / 2 2 . > 

1 -Tr2^_ / w c . . me „ 

-2' ,VSC (^r'B + lV C n iq/ 


(in-18) 


where 


-1 


i B = SCpSc 3 ) B 


B = aircraft moment of inertia about the y-axis. 


The detail derivation of the impedance matrix Z_(w) and the definitions of the 

XV 

aerodynamic derivatives in (H[-18) maybe found in Chapter 6, Reference 55. In the 
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present ease, because of lack of tail surfaces dynamic load data, only the rigid body 
plunging mode is used. In other words, Z R (u) in Equation (HI- 16) is replaced by 

For a free-free mode with index i, the modal influence function is: 

, v _ 1 r {^[i^-^trb]- expfifr-^t ra ] [ exp (i^t) 

( r^ *v ii 2nim t J . .[ 2 . .. * . 2 "I W 1 

11 1 (“-“j) “ L +^^<^-( 1 +^+ 1 /^)^ 


(m-i9) 


where 


m i = < M F >££ 

S t = (g F^£ 

When |> three poles exist in integral (III-19), one of which is located on the 

positive real axis of the -plane; the remaining two are located on the upper half 
of the complex plane. 

For the plunging mode, the modal influence function is: 


\ = 1 1 

r°°\ 

|exp [i(<j - - exp [i(w 

' t ‘'l )t ra]f exp < iw l^ 

/RR 71-pVS J 

-00 (w - <J^) 

\SL7_ilE_c 1 

2V \pSc z a) 

lea. + C i 
1 1 zo 



(m-20) 

The three poles of the integral I r for the rigid body mode are located at the 
origin, the positive real axis and the positive imaginary axis of the complex 
^ -plane. 

In Equations (HI-19, m-20), the corresponding influence functions for modal 

2 

accelerations may be obtained by adding a factor (- — )to the two integrands. These 

G 


functions are denoted by^(t, u> l /(■> 



VII. Acceleration Response of the Aircraft 


n n 


w 2 (t) 


/“£ E lyj r ij( t , ^>j [Ax] r 


r= 1 s= 1 


S (w) 

P P 
r‘s 


[AX] 


Ig * (t»w)J {y T } d 


(DI-21) 


where 


P r P s M 


^ IV"’} L 


T7* T/ \ 

\ <“>j 


S P P (o) 
‘sr 


S (w) 

P P 
. *r s 


r s (~w)T = 

p p 

L ^r^s J 


FS P P M 
*r s 


S (-w) 
P s P r 


S (w) 
P P 




I(t.«) 0 

F 


0 


I* (t,«) 
R 


(IH-22) 


Aw = frequency increment = 2n/AT 

For fp 1 free-free mode, the influence function is obtained by carrying out the 
integration: 


t < t 


ra 


C (t,w) = 


i*-* ■ r 2 2* 

f 2Gm^(l+ e f + i/3j) (ibj-w) -a^ 


2a f w 4 exp (iwt) - P f (t r& ,t,w) 

P f(^ r b* ^“Pf^ra* 


t <t<t , 
ra rb 


t rb <t 


(HI- 2 3) 
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where 




a f - r ~ (* 4 •- + 


f ■ 


1/2 


b, = 


w f g f 


f 1 + c £ + i0 £ 


P f (t 0 ,t,o>) = exp(iwt Q ) {- (-a f + ib f -u) (a f + ib^ 2 exp [i(a £ + ib^ (t-y] 


+ (a f + ib f -u>) ( -a f + ib f ) 2 exp [i(-a f + ib f ) ( t - t Q )]} (ID-24) 

and G = gravitational acceleration. 

For the plunging mode, integration of the modified expression (DI-20) yields: 

Ir R^’ ” (m + p Sc) (w - iy)G 

0 t<t 

ra 

exp(iwt) - yi exp(iwt ra ) exp [-y(t - t ra )] } t ra < t <t rfe 

iy{exp(i< J t rb ) exp [-y(t - t rb )]- exp (iut ra )exp [-y(t - t ra )] } t rb < t 




-C 


za 

c 

/_4m \ 

2V 

Use L za) 


(III— 25) 
(ID-26) 
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The infinite frequency integral in expression w (t) may be truncated as follows: 


2 n n 


w 2 (tk ) 2 H H w 


i-i r=l s=l 




[AX]' 


S (u> .) 

P P i 

. *r s 


[ A x] 


+ fOv -v] c ax ] t [ S p r p s <“i > ] m \€ <*k* -“i>j){ YT } 


where 


(m-27) 


= ijhw - lower limit of integration 


u. - Liu = upper limit of integration 
x 2 2 


The first summation sign on the right hand of equation (III-27) covers the frequency 

range (w. , w- ) with increment hw. In the computation based on dynamic buffet 
X 1 l 2 

pressure of Run 2, Flight 871, the following data are used: 
w f /27r = 4.049, 6.522, 8. 378 Hz 


i j . /2n = 1. 465 Hz 


w. /2rr = 20. 0 Hz 


38 frequency values 


dw/2 n = 0.488 Hz 

4 la = °» *ib = t 2a = 2 * 20 sec ’ Sb = ^a = 4 ' 55 sec ’ ^b = ha. = 6 * 90 sec ’ 


t^ = 9. 10 sec 


A computer program flow chart showing the computation procedure is given in Figure 
90. The numerical results based on the above mentioned flight data are presented in 
Subsection 5.4 where the time varying PSD's of C.G. and wingtip accelerations are 
given in Figures 85-87. 
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FIGURE 90. COMPUTER PROGRAM FLOW CHART 
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